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ABSTRACT 


The  purpose  of  this  study  was  to  devise  workable  procedure 
for  determining  the  geodetic  position  of  an  observation  site  from 
orbital  data  and  photographic  plates  of  artifical  satellites. 

Development  of  procedure  depended  heavily  on  published 
information  from  the  Smithsonian  Institute,  Space  Track,  and  the 
U.  S.  Naval  Observatory. 

Along  with  Space  Track  bulletins,  two  or  more  good  photo* 
graphic  plates  of  a  satellite  are  needed  for  a  solution. 

The  time  lapse  between  plates  should  be  at  a  minimum.  This 
recommendation  stems  from  a  conclusion  that  satellites  now  in 
orbit  are  not  entirely  satisfactory  for  geodetic  purposes  because 
their  orbital  elements  are  changing  too  much  from  day  to  day. 
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INTRODUCTION 


The  uses  for  artificial  earth  satellites  are  increasing  and  expanding  into  moi'e  scien¬ 
tific  fields  than  was  first  expected.  The  differing  sizes  and  shapes  of  the  orbits  as  well 
as  the  varying  instrumentation  provide  data  for  a  wide  range  of  scientific  fields. 

Observations  of  a  satellite  can  help  the  geodesist  to  determine  some  very  basic  and 
essential  data.  A  study  of  the  irregularity  of  the  orbit  can  reveal  data  concerning  the 
earth’s  gravitational  field.  From  investigations  of  a  satellite’s  period  and  acceleration, 
more  precise  data  can  be  collected  from  which  to  determine  the  size  of  the  earth.  The 
rate  of  change  of  the  plane  of  the  orbit  will  help  fix  the  flattening  of  the  earth.  Finally, 
the  satellite  provides  the  geodesist  with  a  direct  method  for  connection  of  datums. 

The  field  of  geodesy  can  use  the  satellites  in  two  basic  ways.  One  method  involves 
observing  the  satellite  simultaneously  Fom  two  or  more  stations.  The  timing  must  be  pre¬ 
cise.  Atmospheric  refraction  ceases  to  be  important  uae  to  the  simultaneity  of  exposure. 

The  other  method  involves  knowing  the  orbit  of  the  satellite  and  then  making  obser¬ 
vations  from  any  station  on  the  earth’s  surface.  The  latter  method  will  be  the  one  treated 
in  this  report.  The  requirement  of  the  satellite  for  the  carrying  out  of  this  method  is  large¬ 
ly  in  the  character  of  the  orbit  rather  than  in  special  instrumentation.  A  steady  circular 
orbit  would  provide  the  necessary  conditions  for  a  satelliit.  to  be  used  in  triangulation  lor 
establishing  the  geodetic  position  of  an  observation  site  or  for  the  connection  of  datums. 

Most  geodetic  needs  depend  on  accurate  prediction  of  the  satellite’s  position;  some¬ 
thing  which  is  not  possible  with  an  orbit  of  large  eccentricity.  The  perturbing  factors  such 
as  the  earth’s  gravity  and  atmospheric  drag  are  more  pronounced  on  a  satellite  the  closer 
it  IS  to  the  earth.  On  the  ether  hand,  a  satellite  which  is  a  great  distance  from  the  earth 
may  be  very  difficult  to  track  either  photographically  or  electronically. 

It  is  highly  unlikely  that  a  circular  oibit  could  oe  obtained.  There  is  no  example  of 
any  orbiting  body  traveling  in  a  circular  path  in  the  known  universe.  But  we  can  expect  to 
attain  an  elliptical  orbit  with  a  very  small  amount  of  flattening  which  would  be  approach¬ 
ing  the  character  of  a  circle. 

The  steadiness  of  the  orbit  is  a  necessity  for  geodetic  purposes  because  in  order 
to  make  use  of  the  satellite  as  a  celestial  trinngulation  point,  its  position  must  be  calcu¬ 
lated  from  known  orbital  data.  If  the  orbit  is  changing  greatly  from  one  revolution  to  the 
next,  the  data  applicable  for  computing  the  position  during  one  revolution  will  not  l)e  true 
for  computations  of  any  position  in  tiie  next  revolution. 

When  the  character  of  the  orbit  more  nearly  satisfies  the  geodetic  requirements,  then 
the  remaining  consideration  in  connection  with  photographic  tracking  is  ’/isibility.  The 
camera  has  to  be  able  to  record  the  image  or  the  path  of  the  satellite  on  the  photographic 
plate.  The  size  and  height  of  the  satellite  as  well  as  characteristics  of  the  orbit  will 
determine  the  visibility  photographically. 

Hspecially  designed  equipment  is  bring  made  for  photographic  satellite  tracking  and 


methods  of  plate  reduction  specifically  adapted  to  each  type  of  equipment  have  been  or  are 
being  worked  out.  The  Markowitz  Dual  Rate  Tracking  camera,  the  Schmidt  camera,  and  the 
fiaker-Nunn  Tracking  camera  are  some  of  the  instruments  which  have  been  designed  esjje- 
cially  for  gathering  satellite  data. 

As  is  true  of  most  scientific  equipment,  the  portability  decreases  with  the  precisicm 
of  the  machinery.  Most  of  the  specially  designed  equipment  cannot  be  considered  portable 
field  equipment  in  the  sense  that  it  can  be  transported  and  assembled  in  a  matter  of  days. 
However,  the  Markowitz  camera  is  considered  feasible  for  observation  statituis  maintained 
for  six  months  to  a  year. 

It  is  the  purpose  of  this  report  to  present  a  method  for  finding  the  coordinates  of  a 
position  on  the  earth  from  data  obtained  from  observations  of  an  artificial  earth  satellite. 
Included  in  the  intention  is  the  recognition  of  the  need  for  a  method  which  is  usable  with 
portable  equipment. 

The  presentation  is  organized  in  three  phases.  The  first  phase  presents  the  formulae 
and  procedure  for  computing  the  coordinates  of  the  satellite  from  published  orbital  data. 
Agencies  with  the  assignment  of  computing  the  orbits  of  satellites  collect  information  from 
their  various  field  stations  and  compute  and  publish  the  necessary  orbital  elements  such 
as  the  inclination,  right  ascension  of  the  ascending  node,  the  distance  at  apogee  and 
perigee,  et  cetera,  from  which  the  coordinates  may  be  obtained.  Smithsonian  Astrophysical 
Observatory  and  Air  Force  Cambridge  Research  Center  are  the  two  agencies  which  supply 
the  orbital  data  and  prediction  information  to  ACIC. 

The  second  phase  describes  the  adaptation  of  classical  methods  of  astronomical 
photography  for  the  purpose  of  obtaining  the  coordinates  of  the  satellite  from  a  photograph* 
ic  plate.  It  is  not  the  intention  of  this  report  to  describe  the  method  of  photographing  the 
satellite  or  recording  the  time  of  observation.  Phase  two  begins  with  the  procedure  for 
identifying  the  reference  stars  after  the  photographic  plate  has  been  processed.  It  must  be 
assumed  that  the  operator  is  trained  in  the  use  of  the  measuring  instrument.  The  reduction 
is  prepared  for  readings  from  a  linear  comparator  type  of  plate  measuring  device;  not  for 
readings  of  angular  measure. 

The  third  phase  combines  the  results  of  the  first  and  second  phases  in  order  to  pro¬ 
duce  the  x,y,z  space  coordinates  of  the  observation  site.  These  can  be  convert  d  to  geo¬ 
centric  latitude  and  longitude.  Then  these  geocentric  coordinates  can  be  converted  into 
geodetic  latitude  and  longitude  by  the  standard  formulae  which  has  been  included  in 
Appendix  F. 

These  three  phases  of  the  report  contain  all  the  necessary  working  formulae,  the 
outlined  procedure,  and  sample  computations.  The  derivation  of  formulae  and  additional 
information  which  is  not  essential  to  the  computation  procedure  is  relegated  to  the  appen¬ 
dixes. 
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THE  DETERMINATION  OF  RIGHT  ASCENSION  AND  DECLINATION 
OF  AN  ARTIFICIAL  EARTH  SATELUTE  AT  A  GIVEN  UNIVERSAL 
TIME  FROM  ORBITAL  DATA 


The  orbit  of  a  satellite  traveling  around  the  earth  is  an  ellipse  with  the  earth  as  one 
of  the  foci.  This  elliptical  path  is  undergoing  continuous  changes  due  to  gravitational  pull 
from  the  earth,  atmospheric  drag,  and  non*sphericity  of  the  earth.  The  effects  of  these 
forces  on  the  satellite’s  orbit  are  referred  to  as  perturbations.  The  result  of  these  com¬ 
bined  perturbing  factors  can  be  determined  by  observing  what  has  actually  happened  to  the 
orbit  of  a  satellite  by  comparing  successive  revolutions.  But  it  is  much  harder  to  predict 
what  their  effect  will  be  on  future  revolutions  of  the  satellite.  Other  elements,  too,  maKe 
predictions  more  difficult;  the  proximity  to  the  earth,  the  ellipticity  of  the  orbit,  and  the 
speed  of  the  satellite.  The  closer  the  satellite  travels  to  the  earth,  the  more  its  orbit  is 
affected  by  the  earth’s  gravitational  attraction;  also  the  lower  it  is  in  the  earth’s  atmos¬ 
phere  the  more  it  is  perturbed  by  the  atmospheric  drag.  A  very  eccentric  orbit  causes  the 
satellite  to  travel  close  to  the  earth  in  part  of  the  orbit  and  many  times  farther  away  from 
the  earth  in  the  opposite  part  of  the  orbit.  The  disturbing  effects  are  not  constant,  then, 
throughout  the  path.  More  disturbances  are  influential  at  the  close-in  point,  called  perigee, 
than  at  the  farthest-away  point,  the  apogee.  The  greater  the  speed  of  any  space  vehicle, 
the  less  its  path  is  diverted  from  stability.  (The  speed  can  vary  from  the  minimum  which 
just  barely  overcomes  the  force  tending  to  draw  a  mass  back  to  earth  to  a  maximum  value 
w!iich  allows  the  space  object  to  completely  overcome  the  earth’s  gravitational  pull,  and 
escape  from  an  orbit  around  the  earth.) 

Many  other  forces  with  lesser  effects  also  contribute  to  the  difficulty  of  determining 
orbital  elements  which  will  hold  true  for  more  than  just  a  few  revolutions.  It  is  necessary 
to  realize  that  the  orbital  elements  are  undergoing  continual  change,  when  the  position  of 
the  satellite  at  a  particular  time  is  to  be  computed  from  the  orbital  data.  If  a  position  of 
the  satellite  is  desired  for  a  particular  time  which  relates  to  the  1500th  revolution,  for 
example,  the  coordinates  of  the  position  will  be  more  dependable  if  the  orbital  elements 
refer  to  that  same  revolution  than  if  the  orbit  was  computed  from  the  1480th  revolution. 

How  much  the  orbit  varies  from  revolution  to  revolution,  of  course  depends  on  the  charac¬ 
teristics  of  the  particular  satellite. 

Care  should  therefore  be  taken  to  select  observations  of  the  satellite  which  have 
been  made  during  (or  as  closely  as  possible  to)  a  revolution  for  which  orbital  data  was 
computed. 

In  order  to  compute  the  coordinates  of  the  satellite  from  orbital  data  we  have  to  con¬ 
sider  the  elements  involved. 


Tim* 

When  the  sun  is  on  the  meridian  of  a  place  it  is  apparent  noon  there;  when  the  sun  is 
next  on  the  meridian,  an  apparent  solar  day  is  said  to  have  elapsed.  The  meridian  of 
Greenwich  is  used  as  the  reference  point  or  the  0  meridian  and  the  solar  time  of  the 
Greenwich  meridian  is  referred  to  as  Greenwich  Civil  Time  (GCT). 
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Starting  from  Greenwich  the  rest  of  the  world  is  divided  into  24  zones.  In  each  zone, 
a  standard  time  appropriate  to  a  definite  meridian  within  the  zone  is  kept.  This  standard 
time  is  also  referred  to  as  zone  time.  This  local  zone  time  can  be  converted  to  GCT  by 
adding  or  subtracting  the  value  of  the  longitude  of  the  central  meridian  in  time.  For  exam¬ 
ple,  the  central  meridian  of  Central  Standard  Time  is  90®W.  which  is  equivalent  to  ♦6*'. 
GCT  is  then  the  local  zone  time  plus  5 'I  (West  longitude  is  added;  East  longitude  is  sub¬ 
tracted) 

By  agreement  GCT  is  also  called  Universal  Time  (UT).  In  referring  to  the  satellite. 
Universal  Time  is  us'jd.  UT  is  counted  from  0  to  24  hours,  beginning  at  Greenwich  mid¬ 
night. 


Observation  date  (or  epoch)  is  usually  expressed  in  decimal  parts  of  a  day.  For  ex¬ 
ample  an  observation  date  of  August  21,  1958  at  10  hours  23  minutes  15.65  seconds  in  the 
morning  would  be  written: 

1958  August  21  10''23"’ 15?65 
or 

1958  August  21.432820 

ElemMts  of  the  Colostiol  Sphoro 

The  reference  surface  or  celestial  sphere  is  considered  to  be  a  sphere  of  infinite 
radius  whose  center  coincides  with  the  center  of  the  earth.  All  astronomical  bodies  are 
thought  of  as  being  located  on  the  celestial  sphere,  regardless  of  their  distances  from  the 
earth’s  center.  It  is  fixed  in  space,  and  since  the  earth  rotates,  the  celestial  sphere  ap 
pears  to  rotate  from  east  to  west. 

The  north  and  south  celestial  poles  represent  the  intersections  of  the  earth’s  axis  of 
rotation  with  the  celestial  sphere.  The  extension  of  the  plane  of  the  earth’s  equator  inter¬ 
sects  the  celestial  sphere  in  the  great  circle  called  the  celestial  equator. 

A  celestial  meridian  is  an  arc  on  the  celestial  sphere  which  is  the  extension  of  any 
terrestial  meridian. 

The  ecliptic  is  the  great  circle  on  the  celestial  sphere  which  is  the  apparent  annual 
path  of  the  sun. 

The  plane  of  the  ecliptic  is  at  an  angle  with  the  plane  of  the  equator.  This  angle  is 
called  the  obliquity  of  the  ecliptic  and  is  approximately  23  1/2  degrees.  The  two  points 
where  the  ecliptic  and  equator  intersect  are  called  the  equinoxes. 

The  vernal  equinox  (T)  is  the  point  at  which  the  path  of  the  sun  crosses  the  equator 
going  from  the  southern  to  the  northern  side. 

To  describe  the  position  of  any  celestial  object  on  the  celestial  sphere  we  use  a 
system  of  coordinates  which  relate  the  position  to  the  vernal  equinox  and  the  celestial 
equator.  Right  ascension  (a)  is  the  angular  distance  measured  eastward  from  the  vernal 
equinox  along  the  celestial  equator  to  the  point  where  the  celestial  meridian  through  the 


4 


object  intersects  the  equator.  The  declination  (8)  is  the  angular  distance  from  the  equator 
to  the  object  measured  along  the  meridian. 

Declination  is  never  greater  than  90**,  but  it  is  positive  north  of  the  equator  and  neg¬ 
ative  south  of  the  equator.  Right  ascension  is  positive  and  is  generally  given  in  terms  of 
time.  The  relationship  between  time  and  arc  is: 

-  15** 
r-  15' 

1*  -  15" 

The  hour  circle  is  a  great  circle  intersecting  the  celestial  equator  at  right  angles  and 
containing  the  north  and  south  celestial  poles  and  the  celestial  object.  It  can  be  consid¬ 
ered  the  same  as  the  celestial  meridian.  It  is  apparent  that  every  celestial  object  will  have 
an  hour  circle. 


Celeatial 
North  Pole 


Fig.  1 

Elements  Pertaining  to  Elliptical  Ori)its 

As  previously  stated,  the  path  of  a  satellite  around  the  earth  is  an  ellipse  with  the 
earth  as  one  of  the  foci.  The  apogee  and  perigee  distances  are  respectively  the  maximum 
and  minimum  distances  of  a  satellite  from  the  earth. 
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Referring  to  Fig.  2:  F  represents  the 
earth’s  center;  A,  the  point  of  perigee;  A', 
the  point  of  apogee.  The  imaginary  line  A 'A 
is  the  major  axis  of  the  orbit.  One  half  of 
A*  A  is  the  semimajor  axis  designated  by  the 
letter  a.  If  S  is  the  position  of  the  satellite 
at  any  point  in  the  orbit,  then  the  distance 
F  S  is  the  radius  vector  R.  Since  the  orbit 
is  elliptical,  the  value  of  R  will  not  be  con¬ 
stant;  therefore,  in  order  to  designate  the 
position  in  the  oiint  to  which  any  particular 
R  refers,  we  must  specify  the  Universal  Fig.  2 

Time  of  the  position.  Two  angles  to  consider  at  this  time  are  the  true  anomaly  (v)  and  the 
mean  anomaly  (M).  Angle  AFS  is  the  true  anomaly  or  the  angular  distance  of  the  satellite 
from  perigee  measured  along  the  orbit.  If  the  angle  AFS  would  be  visualized  as  an  angle 
that  is  zero  at  perigee  and  increases  uniformly  at  a  rate  of  360  degrees  per  orbital  period, 
it  would  then  be  correctly  called  the  mean  anomaly. 

Referring  to  Fig.  3: 

F  is  again  the  earth’s 
center;  S  and  A  are  projections 
of  the  satellite  and  perigee 
point  on  the  celestial  sphere; 

D  is  the  intersection  of  the 
satellite’s  meridian  (hour  circle) 
with  the  celestial  equator;  N  is 
the  projection  on  the  celestial 
equator  of  the  ascending  node 
where  the  satellite,  going  in  a 
northward  direction  in  its  orbit, 
crosses  the  equator. 

Two  elements  specify  the 
orientation  in  space  of  the  orbit¬ 
al  plane;  the  inclination,  i,  and 
the  right  ascension  of  the  ascend¬ 
ing  node  designated  with  the  Greek 
letter,  capital  Omega  (Q).  The 
inclination  is  the  angle  between 
the  orbital  plane  and  the  earth’s  equatorial  plane.  Q  is  the  arc  (T  N),  measured  eastward 
along  the  celestial  equator,  from  the  vernal  equinox  (T)  to  the  ascending  node  (N). 

Another  orbital  element,  ’’the  argument  of  the  perigee”  («),  a  small  omega,  is  the 
angle,  as  seen  from  F,  from  N  to  the  perigee  point.  The  angle  a  lies  in  the  orbital  plane, 
and  is  measured  in  the  direction  of  the  satellite’s  motion. 
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Finally,  three  other  elements,  the  orbital  period  (P),  the  “epoch  of  perigee”  (T), 
and  the  eccentricity  of  orbit  (e)  must  be  mentioned.  If  the  orbit  is  ideal  and  unchanging, 

P  is  simply  the  interval  between  successive  passages  through  any  selected  fixed  point  of 
the  orbit.  In  this  phase,  the  selected  point  will  be  the  ascending  node  and  the  period 
(nodical)  will  be  determined  for  a  specific  time  (UT).  Because  the  orbits  will  not  be  ideal, 
the  change  in  the  nodical  period  with  respect  to  a  change  in  time  will  be  roughly  computed 
representing  the  average  daily  change.  T  is  the  time  of  passage  through  perigee,  e  along 
with  a  are  the  two  elements  which  specify  the  size  and  shape  of  the  orbit  (see  Appendix 
A). 


The  true  anomaly  (v),  in  Hg.  3,  would  be  angle  APS,  pirovided  you  consider  the  angle 
has  not  been  increasing  uniformly.  Arc  SD  represents  the  declination  (fi)  and  the  arc  T  D, 
the  right  ascension  (ac). 

DevelopmMt  of  Formulcw  with  Perturbotiops  Not  Coosidorod 


Once  the  elements  a,  e,  i,  O,,  w,,  T,  and  P  are  known,  the  nt,  Rt,  v, ,  S,  and  a  can 
be  predicted  for  a  particular  time  (t).  n  is  the  mean  angular  velocity  of  a  satellite  in  its 
orbit  expressed  in  radians  per  second.  A  value  for  n  is  given  by: 


"1  = 


Pt 


,  see  Appendix  B,  part  3. 


An  approximate  formula  for  R,  would  be: 


Rt  *  a  f  1-e  cos  n,  (t-T)  1,  as  given  by  Davis,  Whipple,  and  Zirker*. 
T  must  be  the  time  of  the  perigee  passage  nearest  to  t,  time  wise,  but  preceding  t.  Then 
for  V, : 


cos  V,  ■ 


a(l-e")  -  Rt 
e  R, 


the  polar  equation  for  an  ellipse,  (see  Appendix 
A,  part  Id). 


Next  the  declination  and  right  ascension  are  given  by: 

sin  8  -  sin  i  sin  (a>,  -f  v,  ),  where  a>,  and  v,  indicate  values  of  »  and  v  for 
an  observation  time  t. 


a  ^  Q,  +  arc  tan  [cos  i  tan  (ai,  Vt)  ],  see  Appendix  A,  part  2a  and  b  for 
the  basic  spherical  trigonometry  formulae. 


It  should  be  emphasized  that  the  two  formulae  just  given  apply  only  for  an  observer 
at  the  center  of  the  earth,  or  to  one  directly  beneath  the  satellite. 

DtvwlopRWfit  of  Fomwlat  with 
Perturbotions  Considtrtd 

» 

It  is  well  known  that  the  Earth’s  oblateness  and  its  atmosphere  cause  the  main  per¬ 
turbations  to  the  orbit  of  a  near  satellite,  and  various  theoretical  studies  of  these  effects 
have  been  made.  Fortunately,  it  happens  that  the  first  order  perturbations  due  to  the  atmos¬ 
phere  and  the  Earth's  oblateness  are  of  different  kinds;  the  air  drag  steadily  reduces  the 
length  of  the  major  axis  and  the  eccentricity  of  the  orbit,  while  the  Earth’s  oblateness 
causes  the  orbital  plane  to  rotate  about  the  Earth’s  axis  and  the  major  axis  of  the  orbit  to 
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rotate  in  the  orbital  plane.  The  rotation  of  the  major  axis,  expressed  by  changes  in  oi,  is 
difncult  to  determine,  but  the  rotation  of  the  orbital  plane,  expressed  by  changes  inO,  can 
be  measured  very  accurately. 

To  determine  how  much  the  right  ascension  of  the  ascending  node  tf})  changes  in 
radians  each  second  of  time,  the  following  expression  is  taken  from  Merson  and  King*Hele:* 

Change  in  Q  -  n  cos  i  [  J  ♦  sin^  i  - 1)“^  D  B^(7  sin^  i-4) 

14 

In  the  expression  above: 

n , 

a(l-e2) 

where: 


r,*  Earth’s  equatorial  radius 

J,  D  ■  Constants  which  are  functions  of  the  Earth's  oblateness  and  centrif* 
ugal  and  gravitational  acceleration  at  the  equator.  If  based  on  an  oblateness  of 

J  -0.001637 

D  -  10.6  X  10'^,  see  Appendix  C  for  further  explanation  61  J  and  D. 

Information  published  by  the  Astrophysical  Observatory  of  the  Smithsonian  Institu* 
tion  in  Cambridge,  Massachusetts  in  the  form  of  “Announcement  Cards’’  and  bulletins 
published  by  Space  Track  will  list  a  simpliHed  formula  for  determing  the  rate  of  change 
in  for  a  particular  satellite.  The  example  in  this  phase  includes  such  a  formula. 


As  has  been  previously  stated,  when  computing  the  changes  in  u  (argument  of 
perigee)  it  amounts  to  determining  the  extent  of  rotation  of  the  major  axis  of  the  orbit  in 
the  orbital  plane  which  is  beyond  the  scope  of  this  paper. 


The  rate  of  change  in  u  as  found  in  Space  Track  bulletins  incorporates  the  time  of 
the  passage  through  the  ascending  node  (Tf| ).  As  an  example  of  how  this  is  presented,  the 
following  statement  is  taken  from  Space  Track  bulletin  jiQ.58  for  the  satellite  1958  Delta  II: 

“Argument  of  Perigee  287.618  Degrees  Minus  .3668  Paren  T  Minus  Epoch 

Paren.’’ 

To  further  explain  this  statement,  it  should  be  written  in  mathematical  form,  as  follows: 

-  Wo  ~  0.3668(T^  *’  to),  the  T  in  the  statement  is  actually  T^  and  not 
the  perigee  passage,  where 

Wt  -  the  argument  of  perigee  which  is  nearest,  time  wise,  to  an  observa* 
tion  time  t 

Wo  -  287?618,  the  value  computed  for  the  time  to 

Tfj  =  the  time  of  a  passage  through  the  ascending  node  N  which  is  near* 
est,  time  wise,  to  t 

to  =  called  “epoch’’  in  the  statement,  refers  to  that  observation  time  for 
which  the  elements  Qe,  ao,  eo,  Po,  etc.  were  computed 
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The  formula  for  R,  given  by  Moulton’  is: 

R,  -  a(l  —  e  cos  E,),  the  development  of  this  formula  is  given  in  Appen¬ 
dix  B,  part  4. 

where 

E,  -  n,  (t  -  T)  +  et  sin  n,  (t  —  T)+  ^  sin  2  n,  (t  -  T) 
or 

Et  »  M,  +  e,  sin  M,  +^-sin  2  M, 

The  first  term  Mt  must  be  expressed  in  radians.  E  is  called  the  “eccentric"  anomaly  (see 
Appendix  B,  part  4). 

The  formula  for  v,  given  by  Sterne’  is: 

tan  1/2  Vt  tan  1/2  Et 

DtItffniiMtiM  of  Hio  Pofiod 

Space  Track  includes  in  its  satellite  bulletins  a  tabulated  listing  of  the  Greenwich 
meridian  time  (U.T.)  for  the  south-north  equatorial  crossings.  This  listing  may  include  the 
crossings  over  a  period  of  four  to  Hve  days.  Since  these  times  are  for  the  south-north 
crossings  of  a  satellite,  they  refer  to  the  passage  through  the  ascending  node. 

When  the  time  for  a  passage  through  the  ascending  node  is  used  in  a  formula,  it  is 
usually  denoted  by  T^  (which  is  similar  to  the  symbol  T  used  for  the  time  of  a  passage 
through  perigee). 

To  determine  the  period  Pt  of  a  satellite  for  a  particular  observation  time  t,  refer  to 
the  tabulated  listing  and  locate  the  T^  nearest  (time-wise)  to  t.  Calling  this  the  “center" 
T,4,  find  the  difference  between  it  and  the  preceding  T^  and  also  the  difference  between 
it  and  the  next  Tf^.  The  best  value  for  P,  would  then  be  the  average  of  the  two  differences. 
The  period  found  in  this  manner  refers  to  the  nodical  period. 

Procadvrt 


This  section  will  provide  an  outline  and  examples,  based  on  Space  Track  data,  to 
follow  for  determining  the  right  ascension  and  declination  of  an  artificial  satellite,  as  seen 
from  the  center  of  the  earth,  for  a  given  observation  time  t  (U.T.). 

1st  -  List  the  known  orbital  data  such  as  a^,  e^,  i,  and  Ro  which  will  be 
given  for  a  particular  time  t^. 

2nd  -  List  the  formulae  for  determining  at,  e,,  <•>,  and  Rj  foe  time  t.  Also  list 
the  nearest  nodal  passage  (T^^ )  to  time  t  and  list  Pt . 

The  formulae  for  at  and  e,  will  b;  in  the  form  a^  ±  the  daily  change  and 
e^  t  the  daily  change  (see  the  example).  The  daily  change  in  these  elements  is  usually 
minus  and  may  be  determined  by  use  of  two  or  three  successive  bulletins  which  have  these 
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elements  listed  along  with  the  corresponding  times  (tg). 

The  formulae  for  u,  and  12,  involve  and  t^  and  must  be  taken  from  the 
same  bulletin  which  gave  and  and  Qq  being  the  first  term  of  each  formula). 

3rd  -  Tabulate  the  following  terms:  t  -  to,  T^  -  to,  (T,;,  -  to)  ,  a, ,  e, ,  e,^  , 
r  1  +  e,  ,  “ 


1  -  e, 


,  sin  1,  cos 


3  i.  and  n, 


4th  -  Compute  fl, ,  w, ,  and  v, .  (for  v,  use  the  formula  given  by  Sterne) 


The  computation  of  v,  will  involve  computing  values  for  T,  M,  and  E, 
in  that  order.  T,  once  again,  refers  to  the  U.T.  of  the  passage  through  perigee  that  is 
nearest  (time>wise)  to  time  t.  T  is  not  stated  in  Space  Track  bulletins,  but  a  simplified 
method  of  computing  T  is  a  part  of  the  examples  which  follow  this  procedure  outline.  (If 
Smithsonian  bulletins  are  available,  they  will  include  a  direct  formula  for  finding  T.  The 
only  unknown  in  this  particular  formula  is  the  number  of  revolutions  which  the  satellite  has 
made  between  the  time  used  in  the  formula  and  time  t.) 

5th  •  Compute  R,  using  the  formula  given  by  Moulton. 


6th  -  Compute  S  and  a 

sin  8  =  sin  i  sin  (w,  +  v, ) 

a  =  Q,  +  arc  tan  [cos  i  tan  (u^  +  v, )] 
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EXAMPLE 


Given: 


1st  -  ao*’  I.128V20  earth  radii  (this  means  1.128720  times  an  accepted  value  for  the 

equatorial  radius  of  the  earth) 

eo“  0.085822 
i  *  65“  12' 
to  =  Aug.  24.912810  UT 
28?876 
flo*  1057771 
Po*  101.50  minutes 

2nd  -  t  =  Aug.  25.0774535 

*  Aug.  25.0530313  (nearest  nodal  passage  to  t  -  time  wise) 
a,  «  ao  -  0.000442  (t-t^)  (average  daily  change) 
e^  *  Co  -  0.000358  (t-t®)  (average  daily  change) 
u,  *  -  0.350  (T^  -  to)  (from  bulletin) 

fit  “fio  -2.782  (Tn  -  to)  -  0.0018  (T^  -  t,)* 

P,  *  101.45  minutes 

3rd  -Tabulation: 


t  -  to-  0.1646435 
Tn  -  to-  0.1402213 
(T^  -  to)'-  0.019662 
sin  i  =  0.907777 
cos  i  -  0.419452 


a,  -  1.128647  earth  radii 
e,  -  0.085763 
e,  2  -  0.007355 

[4^1**  -  1.089778 
1-e, 


T)  «=  2w/P  (mean  angular  velocity) 

V,  -  6.283185/6.087  x  10^  sec  -  1.032230  x  10'’  rad /sec. 
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4th  -  Computation  of  Q^,  w,,  and  v, 

0,  -  105^771  -  0.390  -  0.000  -  103!381 

*  105*  22'  52"  (see  formula  in  2Dd  step  above) 

-  28f876  -  0.049  -  28?827 

-  28*  49'  3?' 

tan  1/2  V  -  [  -1*  tan  1/2  E, 

*  1-e  ^ 


where; 


Mj  et  sin  M,  -f  et 
%  (t-T) 


sin  2  Mt 


If  "  28*  49’  37",  then  perigee  must  be  28?827  beyond  the  ascending  node  because 
uis  the  angular  distance  of  perigee  from  the  ascending  node  in  the  direction  the  satellite 
is  moving. 


If  *  101.45  minutes,  the  satellite  travels  through  a  mean  angular  distance  of 
3?549/min. 


The  time  (T,^ )  for  the  passage  through  the  ascending  node  is  Aug.  25.0530313;  then 
the  time  T  for  the  passage  through  perigee  is: 


.p  .p  28.827  ..  .  V 

■3Tm9x'60x  24  <■" 


Aug  25.0586719 


Then: 


t-T  =  0.018782  days 

M,  =  1.032230  X  10-5  ^  0.018782  x  24  x  3600 

”  1.675066  radians  or  95*  58'  27"  (mean  anomaly) 

sin  M,  =  0.994569  2  M,  -  191*  56’  54" 

e,  sin  M,  =  0.085297  sin  2  M,  =  -0.207030 

e,  £*-  sin  2  M,  -  -0.000761 

E,  =  1.675066  +  0.085297  -0.000761 

=  1.75%02  radians  or  100*  49'  04"  (eccentric  anomaly) 

1/2  E,  =  50*  24'  32" 
tan  1/2  E,  =  1.209174 
tan  1/2  v,  .  1.317731 
1/2  v,  *  52*  48'  21" 

=  105*  36'  42"  (true  anomaly) 
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5th  -  The  computation  of  R  may  be  included  here  since  it  is  used  in  phase  III 
R,  =  a  (1-e,  cos  E, ) 
cos  E.  =  -0.187686 
e,  cos  E^  =  -0.016097 
1-e,  cos  E,  =  1.016097 

R,  =  1.146815  earth  radii 

6th  -  Computation  of  5  and  a 

sin  5  =  sin  i  sin  (<Jt  v, ) 
w,  +  V,  134°  26'  19" 
sin  (w,  +  V,)  =  0.714001 

sin  8  =  0.648154  (a  negative  sin  8  would  indicate  a  position  below  the 
celestial  equator) 

8  =  40°  24'  09" 

a  =  12,  +  arc  tan  [  cos  i  tan  (w,  +  v, )  ] 
tan  («,  +  V, )  =  -1.019791 
cos  i  tan  (w,  +  v,)  *  -0,427753 

a  *105°  22'  52" +  156°  50’  28" 
a  =  262°  13'  20"  or  28'"  53!33 
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THE  DETERMINATION  OF  RIGHT  ASCENSION  AND  DECLINATION 
OF  AN  ARTIFICIAL  EARTH  SATELLITE  FROM  PHOTOGRAPHIC  OBSERVATION 


A  photographic  observation  of  a  satellite  should  consist  of  a  glass  plate  showing 
the  image  of  the  satellite  against  a  Held  of  stars  and  the  record  of  the  time  corresponding 
to  the  instant  of  exposure.  It  is  preferable  that  the  exposure  show  the  satellite  and  stars 
as  points  and  not  as  trails.  However,  only  the  high  precision  instruments  are  capable  of 
stopping  the  motion  of  the  satellite;  most  field  equipment  is  likely  to  produce  a  plate 
showing  the  stars  as  trails  as  well  as  the  satellite.  It  should  be  remembered  in  evaluation 
of  results  that  point  images  produce  better  accuracy  than  trails  by  a  ratio  of  10  to  1. 

Identification  and  Salection  of  Rtftrtnct  Stors 

At  the  time  of  the  photographic  observation,  some  record  of  the  direction  and  the 
general  position  of  the  path  of  the  satellite  must  be  made  in  order  that  the  reference  stars 
can  be  identiHed.  If  it  can  be  noted  that  the  satellite  was  traveling  through  a  certain  con¬ 
stellation  or  close  to  a  readily  identifiable  group  of  stars  at  the  time  of  exposure,  this 
information  will  be  used  to  identify  the  stars  which  appear  on  the  plate. 

A  more  desirable  method  of  recording  the  portion  of  the  sky  where  the  observation  is 
being  made,  is  to  take  a  series  of  wide  angle  photographs  which  will  include  all  of  the 
exposures  of  the  plates  to  be  used  in  the  reduction. 

If  star  charts  are  available,  they  will  be  very  helpful  in  identifying  the  stars  on  the 
plates.  The  charts  together  with  star  catalogs  can  be  used  to  compare  relative  positions 
and  work  out  the  identifications. 

If  star  charts  are  not  available,  a  small  graph  of  the  area  of  the  photograph  can  be 
made  and  the  positions  of  the  stars  plotted  from  their  coordinates  in  the  catalog. 

A  minimum  of  three  identified  stars  on  the  plate  is  needed  to  obtai..  a  solution  for 
the  position  of  an  unknown  point  on  the  plate,  however,  the  solution  is  improved  by  a 
gn  ater  number  of  reference  stars.  The  recommended  number  of  identified  stars  is  ten.  The 
stars  should  be  selected  so  that  they  are  evenly  distributed  over  the  plate.  They  should 
also  be  selected  so  that  the  image  of  the  satellite  falls  within  the  field  of  the  reference 
stars.  The  solution  is  still  possible  if  the  unknown  falls  outside  the  field,  but  the  results 
may  be  less  accurate. 

The  best  stars  for  the  purpose  are  ones  which  produce  small  and  clearly  defined 
points  of  light.  These  are  easier  to  measure  than  large  stars  which  show  up  as  irregular 
areas  of  varying  density.  The  stars  of  magnitude  1  through  5  are  nearly  always  too  large 
to  measure  accurately.  This  depends  on  the  length  of  exposure  of  the  plate,  but  generally 
the  stars  of  6th  or  7th  magnitudes  provide  the  best  reference  points. 

Roduction  ta  Epoch  of  tha  Co^rdinatts  of  the  Reference  Stors 

From  the  star  catalogs,  find  and  record  the  right  ascension  (o)and  the  declination 
(6)  of  each  reference  star.  These  values  must  be  corrected  for  epoch  and  proper  motion. 
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For  a  star  given  in  the  Boss  catalogs,  epoch  1950,  the  following  formulae  will  give 
the  correction  to  mean  place  for  the  year  of  observation.*^ 


-  a,„  +(10  -  1950)  An.  Var.  ♦  1/2  (to  -  1950)'  1/100  Sec.  Var. 
«o  -  fim  ♦  (to  -  1950)  An.  Var.  +  1/2  (to  -1950)*  1/100  Sec.  Var. 


Where  and  are  the  values  given  for  the  beginning  of  the  epoch,  and  to  is  the  date  of 
the  exposure  in  terms  of  the  year.  For  an  observation  done  after  1950,  the  term  (to  -  1950) 
will  be  positive. 


Sompl*  Compvtotiofi 


Correction  to  mean  place  for  the  year  1958  of  star  number  19320  from  the  Boss  Cata* 
log.  The  Ri^t  Ascension,  1950  is  14**  16'"52*.427.  The  corrections  per  year  are:  Annual 
Variation,  2.4646:  Secular  Variation,  -  0.(K)40;  and  the  third  term,  ♦  0.008. 


a,n  .  14”  16'"52.*427 

(to-1950)  An.  Var.  -  (8)(^2.4646)  +19.7168 

1/2  (to-1950)'  1/100  Sec.  Var.  -  (0.32)(-0.0040)  -00.0013 

3^  -  (0.0000512)(+0.008)  .0000 

-  14"  17"*  12ri425 

The  Declination,  1950  is  +38®  59*  47V29.  The  Corrections  per  year  are  An.  Var. 
-16.549,  Sec.  Var.  +  0.209,  and  3**t +0.024. 


= 


(to-1950)  An.  Var.  =  (8)(-16.549) 
1/2  (to-1950)'  1/100  Sec.  Var.  -  (0.32)( +0.209) 


'  to-1950  T 
100 


3^ 


(0.0000512)(+0.024) 


38°  59*  47?29 
-2'  12'.'392 
+  .067 

+  .000 

38®  57  *  34'.'%5 


Additional  corrections  using  the  Besselian  star  numbers  give  the  coordinates  of  the 
reference  stars  for  the  observation  time.  An  ephemeris  for  the  year  of  the  observation  is 
necessary  in  order  to  obtain  the  constants  for  the  time  of  observation.  The  American 
Ephemeris  and  Nautical  Almanac  for  1958  gives  the  reduction  formulas  on  page  246.  The 
Besselian  numbers,  the  capital  letters  in  the  formulae,  are  obtained  from  the  table  for  the 
day  and  decimal  part  of  the  day  for  the  observation.  Bessel’s  Star  Constants,  the  lower 
case  letters,  are  determined  from  the  following: 

a  *  3!07342  +  1.33613  sin  Oo  tan  8o  a’  -  20'.'()419  cos  Oq 

b  =  cos  Oq  tan  Sg  b*  *  -  sin  «o 
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T  =  the  time  in  units  of  one  year  from  the  beginning  of  the  Besselian  solar  year.  This 
can  be  determined  from  the  table  of  Independent  Star  Numbers. 

gand/<'  =  the  proper  motion  in  right  ascension  and  declination  respectively  obtained  from 
the  star  catalogs. 

aoand8o=  the  values  computed  from  the  formulae  converted  to  the  year  of  observation.  The 
numerical  data  in  both  of  these  reduction  methods  are  for  1958  observations. 
Every  ephemeris  contains  these  conversion  formulae  for  the  particular  year. 

A2  =  terms  of  short  period  in  A. 
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Sample  Computation  by  Independent  Star  Numbers 


14'’  17"’  12!1425 


f 

= 

+  2.499 

f 

= 

+  0.008 

= 

+  0.0002 

^  g  sin  (G  +  aj  tan  8o 

= 

-  0.8612 

^  h  sin  (H  +  aj  sec 

= 

-  0.6935 

14'’  17"’  13?0950 

fio 

= 

38"  57'  34'.‘965 

1 

= 

+  00.0155 

g  cos  (G  +  aj 

-  08.9042 

h  cos  (H  +  aj  sin  &o 

+  10.9322 

i  cos  8o 

+  05.5520 

38"  57'  4215605 

Determination  of  the  Piote  Center 

The  right  ascension  and  declination  of  the  center  of  the  plate  will  be  designated  as 
Ao  and  Dq.  There  are  three  methods  offered  Tor  finding  the  plate  center,  and  the  selection 
can  be  made  on  the  basis  of  which  methal  best  fits  the  particular  conditions. 

Stor  Method 

A  star  near  the  center  of  the  plate  may  be  used,  If  there  is  an  identifiable  star  with 
a  small  and  well  defined  image  which  appears  near  the  center  of  the  plate,  the  right  ascen¬ 
sion  and  declination  of  this  star,  corrected  for  epoch  and  proper  motion  as  explained  in  the 
preceding  section  ran  be  used  for  the  Ao  and  Do  of  the  plate. 

Displacing  the  center  from  the  true  geometrical  center  of  the  plate  by  a  small  amount 
doesn’t  cause  any  ill  effects  in  the  calculation. 

Coiitroid  Method 

This  method  produces  the  center  of  the  field  of  the  reference  stars  which  may  or  may 
not  be  the  geometrical  center  of  the  plate.  After  the  right  ascension  and  declination  for 
each  star  ire  found,  the  center  may  be  calculated  by  taking  the  averages  in  each  direction. 
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This  can  be  written; 


\ 

•  O 


la. 

n 


n 

where  £  denotes  the  summations  for  i  =  1,2,3, . ,n. 


Proportion  Mothod 


Another  way  to  find  the  center  of  the  plate  is  to  pick  the  geometrical  center  G  on  the 
plate,  and  measure  this  location  when  the  reference  stars  are  measured.  By  computing  a 
ratio  of  the  proportional  differences  between  two  reference  stars,  find  the  Ag  and  values 
as  follows.  (See  Fig.  4) 


•  The  values  of  a,  b,  and  c  are  measured  on  the  plate.  (Line  distances  in  millimeters.) 
G  is  the  geometrical  center  of  the  plate.  The  coordinates  of  G  are  given  as  the  differences 
in  right  ascension  and  declination  between  star  number  5  and  star  number  4.  Find  Oq  and 
So  as  follows; 

a-  =  (■ - )  {a.  -  a.)  (ck)  (1/15  sec  estimated) 

a  +  b  ^ 


where 


8.  =  (- 


a  b 


■)ih  “  8.) 


- r  is  the  line  ratio  in  mm.  between  stars  5  and  4. 

a  +  b 

c  is  the  horizontal  distance  from  the  plate  center  (a^,  SJ  to  the  line  connecting 
stars  5  and  4. 


estimated  is  the  value  of  So  as  read  from  the  plate.  (In  practice  a  graph  is 
made,  keyed  to  the  plate  showing  information  given  on  Fig.  4). 

k  is  a  constant  relationship  of  1  revolution  on  the  measuring  engine  or  1  mm. 
on  the  plate  in  seconds  of  arc. 

is  the  difference  in  right  ascension  between  stars  5  and  4. 

85  -  8^  is  the  difference  in  declination  between  stars  5  and  4. 

After  Og  and  8^  values  are  computed,  we  find  A^  and  Dq. 


Ao  =  -  “o 

Do  =  85  -  80 

The  Ag  and  Dq  values  should  be  rounded  off  to  the  nearest  10*  for  \o  ®nd  the  near¬ 
est  1'  for  Dg. 
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Plat*  MaosurMMfits 


A  measuring  engine  of  good  quality  is  necessary  for  obtaining  the  plate  measurements 
of  the  reference  stars,  the  satellite,  and  the  geometrical  center,  when  the  proportional 
method  is  used  for  arriving  at  the  plate  center.  The  measuring  engine  should  have  an  acci> 
dental  error  of  not  greater  than  5  microns  for  each  screw  and  a  systematic  error  of  0.25 
micron. 

The  plate  should  be  oriented  in  the  engine  so  that  the  positive  x  direction  corre> 
sponds  with  decreasing  right  ascension  and  the  positive  y  direction  corresponds  with 
decreasing  declination.  In  other  words,  the  star  with  the  largest  right  ascension,  should 
have  the  smallest  millimeter  measurement  for  x,  and  the  star  farthest  north  or  having  the 
largest  angle  of  declination  should  have  the  smallest  y  reading. 

The  measurements  for  the  center  of  the  plate  depend  upon  which  system  is  selected. 

If  a  star  is  selected  to  be  substituted  for  the  center,  the  measurements  for  its  coordinates 
must  be  found  along  with  those  of  all  the  reference  stars.  When  the  centroid  method  is  used 
there  will  not  be  any  measurements  for  the  center.  When  the  proportional  method  is 
selected,  the  measurements  must  be  obtained  as  described  in  the  p^-evious  section. 

To  obtain  the  x  and  y  coordinates,  the  center  reading  must  be  subtracted  from  each 
corresponding  reading  of  the  stars  and  satellite.  Where  Xn,  Yn  represent  the  readings  of 
the  stars  and  Xq,  Yq  represent  the  readings  of  the  center  by  any  method,  we  will  find  the 
coordinates  Xn,  y^  by  the  following  relationship: 


When  the  centroid  method  is  used  to  compute  the  coordinates  of  the  center,  the 
notation  for  the  center  is  x^,  yo.  This  is  essentially  translating  the  origin  of  the  x,  y  axes 
to  (0,0). 


Computation  of  th*  Standard  Coordinotos  of  th* 

Roftrenc*  Stcrs 

In  finding  the  relationship  between  the  right  ascension  and  declination  of  stars  and 
their  rectangular  coordinates  on  a  photographic  plate,  it  is  necessary  to  employ  the  use  of 
a  set  of  coordinates  which  provide  the  transition  between  the  celestial  coordinates  and 
the  coordinates  from  the  plate  measurements.  These  transition  coordinates  actually  lie  in 
a  plane  which  is  tangent  to  the  celestial  sphere  and  is  parallel  to  the  camera  plate.  The 
geometrical  relations  and  derivations  will  be  given  in  Appendix  D,  however,  here  it  is 
only  necessary  to  know  that  these  coordinates  are  referred  to  as  the  standard  coordinates, 
i  and  1). 

The  standard  coordinates  are  computed  for  each  reference  star  by  using  the  following 
formulae:* 
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^  _ cot  6  sin  (g  -  Ap) _ 

sin  Do  +  cos  Do  cot  5  cos  (a  -  k^) 

^  ^  cos  Dp  -  cot  8  sin  Dp  cos  (a  -  Ap) 

sin  Dp  +  cos  Dp  cot  8  cos  (a  -  Ao) 

These  values  can  be  checked  by  substituting  £  and  t)  of  each  star  in  the  equation 
which  will  be  presented  in  the  final  sections. 

Computation  of  tho  Standord  Coordinotas 
of  tho  Unknown  Point 

The  next  step  is  to  find  the  standard  coordinates  of  the  unknown  point.  If  an  obser¬ 
vation  could  be  completely  free  from  error,  the  plate  coordinates  would  be  identical  with 
the  standard  coordinates.  Since  an  observation  is  not  without  error,  the  discrepancies  must 
be  considered  in  order  to  make  use  of  the  observational  data.  The  errors  inherent  in  the  x 
and  y  coordinates  are  due  to  refraction,  aberration,  and  mechanical  errors  caused  by  the 
camera  and  measuring  engine.  (A  more  specific  discussion  of  the  errors  appears  in  Appen¬ 
dix  D.)  It  has  been  found  that  these  errors  have  a  linear  effect  on  the  plate  coordinates,  and 
the  relationship  between  the  plate  coordinates  and  the  standard  coordinates  of  any  point 
can  be  expressed:^ 

4-x  =  ax  +  by  +  c 

i}-.y=:dx  +  ey  +  f 

In  these  equations,  a,b,c,d,e,f,  are  small  quantities  which  depend  on  the  errors  in¬ 
volved.  They  are  called  the  plate  constants. 

In  general,  x  differs  from  ^  and  y  differs  from  tj  by  small  quantities,  and  we  can 
write  these  equations  in  an  alternate  form  without  loss  of  accuracy. 

^-x=^a  +  T|b  +  c 
Tj“y=^d  +  Tie+f 

Theoretically,  the  plate  constants  can  be  computed  from  either  set  of  equations,  but 
since  the  {  and  tj  values  can  be  computed  more  accurately  than  the  x  and  y  values  can  be 
measured,  it  is  recommended  that  equations  using  the  standard  coordinates  instead  of  the 
plate  coordinates  be  used  to  determine  the  plate  constants. 

Two  methods  for  solution  are  included  here. 

The  first  method  is  to  obtain  the  specific  values  of  the  six  plate  constants.  We  can 
write  three  equations  with  three  unknowns  which  can  be  solved  by  determinants  or  by  any 
method  for  solving  simultaneous  equations. 

Using  the  form  of  equation  with  the  standard  coordinates  we  can  write  n  equations, 
where  n  is  the  number  of  reference  stars. 
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♦  bi|,  +  c 
^2  -  *  a^2  +  btjj  ♦  c 


“  *n  *  ■*■  b^„  ♦  C 

These  are  called  condition  equations  from  which  we  can  form  the  normal  equations. 
The  procedure  is  as  follows: 

First  we  multiply  each  i  and  t|  by  the  focal  length  of  the  camera  in  millimeters.  It 
was  previously  stated  that  i  differs  from  x  and  from  y  by  small  amounts.  When  i  and  ti 
are  multiplied  by  the  focal  length,  it  is  noticeable  that  the  differences  are  small.  We  then 
set  up  a  schematic  arrangement  as  shown  below: 


where 

Sj  =  f  i  +  1  i  ♦  1 

and 

mi  =  -  X, 

coefficients 

of  a 

of  b 

of  c 

s 

m 

1 

Si 

m, 

• 

’la 

• 

1 

Sa 

• 

m  2 

• 

• 

^  n 

1 

Sn 

nin 

First,  take  each  ^  and  multiply  it  by  each  quantity  in 
column. 

its  row.  Then  find  the 

sum  of  each 

^2^2 

^a’la 

(2 

^  2  ^  2 

^  n 

’In 

^nSn 

u,] 

U.s  ] 

[i.m,] 

The  square  brackets  denote  summations  for  i  =  1,  2,  3,  ...,  n. 
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In  the  same  way,  proceed  with  ij  to  produce 


,  [ti,1  ,  [njSil  ,  and 

Also  find 

[sj  mil  and  [m,] 

If  these  computations  are  correct,  we  should  have 


U|S,]  -  u,£,l  . 

Ujiiil  + 

I^.l 

hi  S,]  = 

Ui^iil  + 

+ 

M 

[s.l  = 

1^1,1  + 

[mil 

[S|m,l  =  U  m,l  + 

h|Hiil  + 

Im.l 

If  these  checks  all  work  out  satisfactorily,  then  we  can  go  ahead  and  form  the  normal 
equations  as  follows; 

+  U.Tj.lb  +  [Uc  = 

+  [T|,T|,lb  +  [iijjc  =  -  Xj)] 

liil  a  +  hil  b  +  nc  =  [(Ci  -  Xj)] 

These  equations  can  now  be  solved  for  a,  b,  and  c. 

The  procedure  for  finding  d,  e,  and  f  is  the  same.  The  only  difference  in  setting  up 
the  schematic  arrangement  is  that 

m,  =  -  y, 

The  resulting  normal  equations  in  d,  e,  and  f  are: 

+  UjTijle  +  [^Jf  =  [^j(Tjj-yj)J 

UjT|j]d  +  [TJjTijle  +  (rjilf  =  hjfrjj  -  yj )] 

Ujl  d  +  hjl  e  +  nf  *  Kn, -Yi)! 

To  check  the  computations  of  the  plate  constants,  substitute  the  values  back  into 
condition  equations. 

After  the  values  of  the  plate  constants  have  been  found,  the  standard  coordinates 
and  T),  of  the  unknown  point  can  be  computed  by  substituting  in  either  pair  of  equations 
which  were  used  to  find  the  plate  constants. 

For  the  same  reason  it  was  recommended  to  use  the  standard  coordinates  for  finding 
the  plate  constants,  it  is  recommended  that  the  following  equations  be  used  to  find  the 
standard  coordinates  of  the  unknown  point. 


-  X,  *  a{,  +  brj,  +  c 
Yt  =  +  eri.  +  f 
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in  these  equations,  since  both  unknowns  appear  in  each  equation,  we  must  solve  them 
as  simultaneous  equations. 

((  and  If,  must  now  be  divided  by  the  focal  length  for  use  in  the  final  equations 
which  are  presented  in  the  section  entitled  Celestial  Coordinates  of  the  Satellite. 

Another  calculation  for  obtaining  the  standard  coordinates  is  called  the  method 
of  dependences.”  Instead  of  solving  for  the  specific  values  of  a,  b,  ...,  f,  we  compute 
certain  quantities  which  depend  on  the  reference  stars  selected  and  on  one  position  of  the 
satellite.  The  development  of  this  method  will  be  more  fully  discussed  in  Appendix  D. 

The  terms  a,  b . f  are  replaced  by  P,  0.  and  R.  The  resulting  equations  which  we 

can  now  solve  are: 

P  [XiXj]  +  Q  Ix^y,]  +  R  [xj  =  X, 

P  [xjyj  +  0  [y.yj  +  R  Cy.l  =  yt 

P[xJ  +  Oty.l  +  Rn  =1 

When  the  centroid  method  of  obtaining  the  plate  center  is  used,  these  three  equations 
become: 


P  [x,Xil  +  0  [xjyil  =  xt 
P  [x,y,l  +  Q  [y,yil  =  y, 

Rn  =  1 

The  fact  that  the  sum  of  the  x  coordinates  and  the  sum  of  the  y  coordinates  are  zero 
provides  a  useful  check  on  the  x  and  y  coordinates  which  are  computed  when  the  plate 
center  is  the  centroid. 

After  P,  Q,  and  R  are  determined,  the  dependences  of  the  reference  stars  are  formed 
by  means  of 

n,  =  Px,  +  Qyj  +  R 

A  check  on  this  calculation  is 

n 

2  D,  -  1 
i  =  1 
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Now  that  the  dependences  have  been  found  we  can  compute  the  standard  coordinates 
of  the  unknown  point. 

For 

n 

-  X,  U,(f, -X,) 

i-  1 

and  for  ij, 

n 

i|,  -  y,  .  2  U.  (i|i-y,) 

i  =  1 

In  this  computation,  just  as  in  the  first  solution,  the  focal  length  must  be  considered. 
The  $  and  tj  coordinates  of  the  reference  stars  can  be  multiplied  by  the  focal  length  (in 
millimeters)  or  the  x  and  y  coordinates  of  the  reference  stars  and  the  satellite  can  be 
divided  by  the  focal  length.  If  the  second  alternative  is  chosen,  and  tj,  will  be  in  the 
proper  units  to  substitute  directly  in  the  formulae  for  finding  the  standard  coordinates 
of  the  unknown  point.  If  the  first  alternative  is  selecteu,  and  ij,  must  be  divided  by  the 
focal  length  before  using  them  in  the  formulae. 

This  method  of  dependences  is  best  suited  to  the  type  of  observation  where  a  series 
of  plates  is  taken  of  the  same  portion  of  the  sky;  i.e.,  that  the  plates  all  have  the  same 
center,  and  the  position  of  the  unknown  point  changes  very  little  between  plates. 
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CtUsHal  Coordinates  of  iho  Sotellite 


The  standard  coordinates  and  ij|  determined  from  either  of  the  two  rethods  are 
now  used  to  find  the  right  ascension  (ni)  declination  iSt)  of  the  satellite. 

First  we  find  oj 

tan  (n,  -  Ao)  =  Y - iT 

1-7?,  tan  Do 

Then  using  the  values  of  a,  from  this  equation,  we  can  find  8,  from: 

_ 4  c  _ /  *  \  ^ 


cot  8,  cos  (n,  -  Ao) 


1?,  +  tan  Do 


c  .  /  t  \  it  sec  Do 

cot  8,  sin  (n,  -  Ao)  =  — - 

I?,  +  tan  Do 

Procoduro 

The  given  information  for  the  computation  is  the  right  ascension  (a)  and  declination 
(8)  of  each  reference  star,  the  plate  measurements  of  the  reference  stars  and  the  unknown 
point,  and  the  focal  length  of  the  camera. 

The  computation  procedure  is  as  follows: 

After  the  identification  and  selection  of  the  reference  stars,  their  coordinates  must 
be  corrected  to  the  ..me  of  observation.  This  can  be  done  by  using  either  the  Besselian 
Star  Constants  or  the  Independent  Star  Numbers.  The  choice  of  the  latter  method  is  recom¬ 
mended  in  the  reduction  of  satellite  plates.  If  there  were  many  observations  on  different 
dates  of  the  same  group  of  stars,  the  use  of  the  Besselian  Star  Constants  would  have  the 
advantage,  but  usually  each  satellite  plate  includes  a  different  group  of  reference  stars, 
and  so  the  Independent  Star  Numbers  method  is  better  suited  to  the  computation. 

The  second  step  is  to  find  the  astronomical  coordinates  (Aq,  Do)  of  the  center  of  the 
photographic  plate  by  the  Centroid  Method.  The  equations  are: 


Third,  find  the  plate  coordinates  (xo  yo)  of  the  center  of  the  plate  by  the  Centroid 
Method  also.  For  the  plate  coordinates  this  can  be  written 


The  capital  letters  X  and  Y  refer  to  the  readings  just  as  they  are  taken  from  the  plate.  The 
lower  case  letters  x  and  y  are  the  plate  coordinates  which  have  been  referred  to  the  plate 
center.  Do  not  include  the  measurements  of  the  unknown  point  X,  and  Y,  in  this  step. 

The  fourth  step  is  to  find  the  plate  coordinates  of  the  reference  stars.  Using  the 
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plate  center  found  in  step  three,  substract  the  x  coordinate  from  each  X  reading;  and  the 
y  coordinate  of  the  center  from  each  Y  reading. 


This  step  can  be  checked  by  the  relationship 
2:x„  .  5:y„-0 

Find  the  plate  coordinates  of  the  unknown  point 
X,  -  X,  -  Xo 
Yt  -  Y,  -  y. 

The  fifth  step  is  to  compute  the  standard  coordinates  ^  and  tj.  Note  that  the  denomi¬ 
nators  are  the  same  in  the  two  equations.  Declination  is  given  in  degrees,  minutes,  and 
seconds  of  arc.  Right  Ascension  is  gi'  en  in  hours,  minutes  and  seconds  of  time,  and  must 
be  converted  to  arc  by  using  the  following  relationship: 

!*•  .  15“^  i*"  -  15'^  1*  .  15" 

After  the  standard  coordinates  have  been  computed,  multiply  each  {  and  i)  by  the  focal 
length  in  millimeters. 

Sixth,  find  the  plate  constants  a,  b,  c,  d,  e,  f.  It  is  recommended  that  the  standard 
coordinates  be  used  for  this  solution  rather  than  the  plate  coordinates. 

|-x»^a+Tjb  +  c 

TJ  -  y  -  4d  +  »je  +  f 

The  solution  for  a,  b,  and  c  is  carried  out  by  following  the  least  square  solution  described 
in  the  section  entitled  ‘‘Computation  of  the  Standard  Coordinates  of  the  Unknown  Point". 
Then  the  equations  for  d,  e,  and  f  are  solved  by  repeating  the  procedure. 

The  standard  coordinates  for  the  unknown  point  are  then  found  by  making  use  of  the 
plate  constants  from  step  six. 

—  X,  »  a  +  TJ,  b  +  c 

»/t  -  Yt  -  d  +  TJ,  e  +  f 

These  equations  have  to  be  solved  simultaneously  for  and  tj,.  When  and  tj,  are 
divided  by  the  focal  length  they  can  be  substituted  in  the  same  equations  used  in  step 
five  to  find  the  right  ascension  and  declination  of  the  unknown  point. 


Somple  Computation 


Coordinates^  corrected  to  epoch,  of  the  reference  stars 
Star  Boss  Catalog  fiumber  a  _ 8 


' 

1 

19429 

u'* 

22*" 

02?2994 

+  41°  02' 

09'.'1446 

2 

19320 

14 

17 

13.0946 

+  38  57 

42.5604 

3 

19225 

14 

12 

42.7228 

+  41  43 

01.8724 

4 

19124 

14 

08 

17.2595 

+  40  58 

38.4310 

5 

19414 

14 

21 

49.9457 

+  39  31 

17.4283 

6 

19322 

14 

17 

17.4187 

+  42  12 

01.9742 

Plate  Measurements 

Star 

X 

Y 

1 

60.4910 

45.7711 

2 

69.6317 

54.0542 

3 

67.6740 

38.4797 

4 

73.4747 

40.1880 

5 

63.9537 

53.2307 

6 

62.3659 

38.0226 

Satellite 

Xt 

Yt 

68.6078 

48.6072 

Focal  length  of  camera 

311.66  mm 

Astro  coordinates  of  the  center  of  the  field  of  stars 


,  85"  39'"  22!7407 

Ao  - - 2 - =  14  16  33.7901 


,,  244“  24'  51V4109  ..o 

Do  * - 7 - =  40  44  8.5684 


Plate  coordinates  of  the  center  of  the  field  of  stars. 
397.5910 


6 

269.7463 
yo — z — 


66.2652 
=  44.9577 
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Plate  coordinates  of  the  reference  stars. 


Star 


Xn=  Xn  -  Xo 


1 

2 

3 

4 

5 

6 


-  5.7742 
+  3.3665 
+  1.4088 
+  7.2095 
-2.3115 

-  3.8993 


+  0.8134 
+  9.0965 

-  6.4780 

-  4.7697 
+  8.2730 
-6.9351 


Plate  coordinates  of  the  satellite 

X,  •=  X,  -  x„ 

’  +2.3426 


Yt  *  Yt  -  Yo 

+3.6495 


Computing  the  Standard  (Coordinates  of  each  reference  star. 

cot  S  sin  (g  -  Ap)  cos  (a  -  Ao) 


+  1.148914781 
+ 1.236580996 
+  1.121697783 
+  1.151287611 
+  1.21216%28 
+  1.102825132 


+  0.023887597 
+  0.002858299 

-  0.016802899 

-  0.036100879 
+  0.022989457 
+  0.003172755 


+  0.999714651 
+  0.999995915 
+  0.999858821 
+  0.999348151 
+  0.999735708 
+  0.999994%7 


sin  Dp 

+  0.652570836 

(311.66)(0 
+  5.616602180 
+  0.693002317 
-  3.909829081 
-8.497551116 
+  5.528992968 
+  0.732758139 


cos  Dp 

+  0.757727725 

(311.66)(q) 

+ 1.676782734 

-  9.651529777 
+  5.361103341 
+  1.414961445 

-  6.565193637 
+  7.970477097 


Condition  equations  for  finding  the  plate  constants  a,  b,  and  c: 

5.616602180  a  +  1.676782734  b  +  c  -  11.390762180 

0.693002317  a  -  9.651529777  b  +  c  .  -  2.673527683 

-  3.909829081  a  +  5.361103341  b  +  c  -  -  5.318659081 

-  8.497551116  a  +  1.414%1445  b  +  c  .  -  15.707081116 

5.528992%8  a  -  6.565193637  b  +  c  -  7.840452%8 

0.732758139  a  +  7.970477097  b  +  c  -  4.632018139 
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-  150.628308411 
[f,  -  X,]  -  0.163%5407 

-  xj]  -  -20.390027237 
IfiV,]  “  -  60.713893397 

233.337444947 

-  y.^l  —  178.104463367 


hi(n,  -  y,)l- 470.836618882 
[{,]-  0.163975407 

-  X,)  1-  263.135349755 
0.206601203 
-  y,  1  -  0.206591203 


Normal  equations  for  finding  a,  b,  and  c: 

150.628308411  a  -  60.713893397  b  +  0.163975407  c  -  263.135349755 

-  60.713893397  a  +  233.337444947  b  +  0.206601203  c  =  -20.390027237 
0.163975407  a  +  0.206601203  b  +  6  c  -  0.163%5407 

Solving  these  equations  by  determinants: 

a-  1.912312557 
b-  0.410229843 
c  -  -  0.039060134 

Condition  equations  for  finding  the  plate  constants  d,  e,  and  f: 

5.616602180  d  +  1.676782734  e  +  f  -  0.863382734 
0.693002317  d  -  9.651529777  e  +  f  -  -18.748029777 

-  3.909829081  d  +  5.361103341  e  +  f  - 11.839103341 

-  8.497551116  d  +  1.414%1445  e  +  f  -  6.184661445 
5.528992968  H  -  6  565193637  e  +  f  -  - 14.838193637 
0.732758139  d  +  7.970477097  e  *■  f  -14.905577097 

Normal  equations  for  finding  d,  e,  and  f: 

150.628308411  d  -  60.713893397  e  +  0.163975407  f  =  -178.104463367 

-  60.713893397  d  +  233.337444947  e  +  0.206601203  f  =  470.836618882 
0.163975407  d  +  0.206601203  e  +  6  f  =  0.206591203 

Solving  by  determinants: 

d  =  -0.412290640 
e  =  1.910576330 
f  =-0.020088440 

Standard  coordinates  of  the  satellite: 

1.912312557fj  +  0.410229843,,,  -  0.039060134  =  -  2.3426 

-  0.412290640  +  1.910576330,,,-  0.020088440  =  ,,,  -  3.6495 

4,  =  -  0.608774171 
q,  =  -  4.261458733 
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Right  ascension  and  declination  of  the  satellite: 

,  .  ,  (-0.608774171  )(1.319735265) 

an  la,  -  A  J  ^  ^  {i.261458733)(0.861220745) 

- -0.002547872 

a,  -  H*'  15'"58!754 

cot  6.  sin  (a,  -  A  )  = -0.(X)3041571 

I  l  O 

S,=  39®  57'  08!«689 
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THE  DETERMINATION  OF  THE  GEOCENTRIC  COORDINATES  OF  A  POSITION  ON 
THE  EARTH’S  SURFACE  FROM  SATELLITE  OBSERVATIONS 


In  the  first  phase  of  this  report  the  {Mxx:edure  is  presented  for  finding  the  coordinates 
of  a  satellite  by  computation  using  the  elements  of  the  orbit.  The  resulting  coordinates 
have  as  their  origin  the  center  of  the  earth. 

In  phase  two  the  method  is  given  for  finding  the  coordinates  of  a  satellite  from  a 
photographic  observation.  The  coordinates  determined  from  observations  are  based  on  an 
origin  located  at  the  observation  station. 

The  coordinates  of  the  satellite  which  are  determined  from  orbital  data  are  those 
which  would  be  seen  by  an  observer  at  the  center  of  the  earth.  These  are  called  geocentric 
coordinates  of  the  satellite.  Those  which  are  observed  from  the  surface  of  the  earth  are 
called  the  observed  or  apparent  positions. 

The  parallax  of  the  satellite  is  the  difference  of  directions  of  the  straight  lines 
drawn  to  the  satellite  from  the  two  different  reference  points.  It  may  also  be  defined  as 
the  angle  at  the  satellite  made  by  the  line  connecting  the  center  of  the  earth  and  the  satel¬ 
lite  and  the  line  connecting  the  observation  station  on  the  surface  of  the  earth  and  the 
satellite.*'^ 

The  observer  at  0  in  figure  5  views  the  satellite  S  along  the  line  OS.  When  the 
center  of  the  earth  is  the  reference  point,  the  position  of  the  satellite  is  referred  to  line 
CS.  The  angle  CSO  is  called  the  parallax  of  the  satellite. 


Fig.  5 

The  coordinates  obtained  from  the  photographic  observation  differ  from  the  ones 
resulting  from  use  of  the  orbital  data.  The  discrepancy  describes  the  difference  in  per¬ 
spective  between  the  earth  centered  location  and  the  position  on  the  earth’s  surface.  If 
the  satellite  were  precisely  at  the  observer’s  zenith,  there  should  be  no  discrepancy 
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between  the  observed  coordinates  and  the  computed  geocentric  coordinates,  and  it  would 
not  be  possible  to  compute  the  observer's  position  from  satellite  observations. 

The  location  of  the  satellite  for  the  photographic  observation  Is  of  considerable 
importance.  The  ideal  location  would  be  a  small  angular  distance  from  the  observer’s 
zenith.  Very  incongruous  results  are  obtained  if  a  satellite  is  photographed  low  on  the 
horizon. 

Since  the  satellite  is  a  great  deal  closer  to  the  earth  than  any  of  the  reference  stars, 
the  projection  of  the  satellite’s  image  onto  the  field  of  stars  depends  on  the  origin  of  th« 
projection.  Suppose  the  satellite  is  in  position  S  in  figure  6,  the  observer  at  K  sees  the 
satellite  projected  on  the  celestial  sphere  at  the  point  marked  observed.  When  the  satellite 
isn’t  far  from  the  observer’s  zenith  the  observed  and  computed  positions  are  close  together. 
The  lower  a  satellite  is  toward  the  horizon,  the  more  its  observed  position  differs  from  the 
computed  position  which  is  referenced  to  the  center  of  the  earth. 

It  is  actually  the  fact  that  there  exists  a  difference  between  the  earth-centered 
position  and  the  observed  position  which  allows  us  to  compute  the  observer’s  position. 
However,  the  formulae  depend  on  the  parallax  angle  being  small. 

The  present  orbits  of  the  satellites  are  undergoing  considerable  changes  due  to  the 
perturbations  already  discussed.  These  changes  in  the  orbit  cause  inaccuracies  in  the 
results  of  phase  one  which  in  turn  affects  phase  two.  Phase  one  is  hindered  if  the  times 
of  the  observations  to  be  used  are  too  widely  separated  or  if  there  is  much  of  a  time  lapse 
between  the  epoch  of  the  furnished  orbital  information  and  the  time  of  the  observation.  It 
is  not  recommended  that  observations  be  any  farther  apart  than  two  consecutive  orbits. 

The  accuracy  loss  which  results  from  the  time  lapse  between  the  furnished  orbital 
data  and  the  time  of  observation  depends  on  each  particular  satellite.  The  orbital  informa¬ 
tion  for  a  satellite  with  a  fairly  steady  orbit  suffers  less  accuracy  loss  between  revolutions 
than  the  satellite  which  is  experiencing  severe  orbital  disturbances. 

Many  other  factors  contribute  to  the  loss  of  accuracy,  but  the  need  here  is  more  to 
point  out  the  fact  that  many  different  factors  can  contribute  to  errors  in  the  results  than  to 
itemize  the  possible  trouble  spots.  Careful  selection  of  the  input  data  for  phase  three  is 
very  essential  to  reasonable  results. 

Methods  for  Finding  the  Observer’s  Cooidinotes 

Several  methods  exist  for  finding  the  position  of  an  observation  station  by  using 
some  celestial  object  as  a  triangulation  point. 

The  Moon  Camera  Method  for  finding  the  observer’s  position  developed  by  Dr.  Wm. 
Markowitz  of  the  L.  S.  Naval  Ob.e.:;  vatory  is  adaptable  for  use  with  an  artificial  earth 
satellite.  I'he  two  methods  which  follow  have  been  adapted  to  fit  the  situation  which 
exists  in  the  case  of  the  satellite.'^  The  Linear  Method  provides  a  way  to  find  the  coordi 
nates  of  an  unknown  station  while  the  Differential  Method  results  in  corrections  to  an 
approximate  |, osition. 
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Linear  Method 


Each  observation  of  the  satellite  yields  the  coordinates  a\  S‘  and  the  time  from  the 
photographic  plate  reduction  phase;  the  coordinates  a, 8  and  the  value  of  K  from  cornpu* 
tation  based  on  orbital  data.  These  six  values  provide  the  input  for  one  observation. 

At  least  two  observations  are  necessary  in  order  to  solve  for  the  geocentric  coordinates 
^01  yoi  of  the  observation  station.  The  equations  which  relate  all  the  observed  and 
computed  data  are 

sin  G'  +  yo  cos  G'  *  R  cos  8  sin  («'  -  a) 
cos  G'  -  Yo  sin  G'  -  z„  cot  5  ‘  “  R  cos  8  cos  (a'  -  a) 

*  R  sin  8  cot 

where: 

G'-  o'-  y 

y  -  The  hour  angle  of  the  Vernal  Equinox  with  respect  to  Greenwich 

R  >  The  distance  from  the  center  of  the  earth  to  the  satellite  in  units  of  earth 
radii. 

In  the  numerical  example  included  in  this  phase,  there  is  an  explanation  of  the 
treatments  of  the  time  of  the  observation  including  the  determination  of  y. 


Since  R  must  be  used  in  the  equations  in  units  of  earth  radii,  the  values  of  Xo,  yb 
and  Zo  will  be  in  the  same  units.  To  Rnd  their  equivalent  values  in  meters,  multiply  each 
by  the  semimajor  axis  (a)  of  the  reference  ellipsoid  of  the  area.  The  numerical  example 
which  follows  gives  us  the  geocentric  coordinates  of  a  station  in  North  America,  and 
therefore,  we  multiply  each  coordinate  by  the  a  of  the  Clarke  1866  ellipsoid. 


ELLIPSOID 

Clarke  18^ 
International 
Hough 
Bessel 


SEMIMAJOR  AXIS -a 

6,378,206.4  Meters 
6,378,388  Meters 

6,378,270  Meters 

6,377,397.155  Meters 


Differential  Method 


Again  we  need  at  least  two  observations  to  get  a  solution.  The  results  to  be  obtained 
are  small  corrections  to  be  applied  to  the  space  coordinates.  First,  we  start  with  an 
approximate  position.  If  the  geodetic  coordinates  are  the  given  information,  the  x,  y,  z 
space  coordinates  need  to  be  computed;  Appendix  H  gives  the  necessary  information  for 
that  conversion.  We  denote  these  preliminary  coordinates  as  x,,  y,,  z,,  and  using  them  in 
the  following  equations,  we  find  and  8^. 
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j  _  X,  sin  G  +  y#  cos  G 

-X,  cos  G  +  Ya  sin  G  +  li  cos  fi 
cos  («c  -  (»)  (|{  sin  8  *  Za) 

X,  cos  G  +  Ya  sin  G  +  H  cos  ^ 

These  values  just  obtained  and  the  observed  values  are  combined  to  find  *\i»' and 
SS\  Both  of  these  increments  should  be  small. 

rt'  -  ftj  =  Aa' 

S'  ~  8'c  =  AS' 

As  long  as  the  increment  of  any  of  the  coordinates  is  small,  it  can  be  considered  the 
same  as  the  differential.  Therefore,  we  can  relate  all  the  observed  data  with  the  unknowns 
by  two  differential  equations  to  be  solved  for  Ax^,,  Ay^,,  AZq. 

a^dXg  +  a^dyo  =  a3da' 

b,dXo  +  badyo  +  b3dzo  =  bade'  +  bjd^' 

where; 

a,  =  sin  G  +  tan  (a'  -  a)  cos  G 
aj  =  cos  G  -  tan  (a'  -  a)  sin  G 

aj  *  sec^  (o'  -  a)  (R  cos  8  -  x,  cos  G  -  Ax  cos  G  +  y,  sin  G  +  Ay  sin  G) 

bj  =  tan  s'  cos  G 

bj  =  -  tan  8  sin  G 

bj  =  -  cos  (o'  -  a) 

b4  =  sin  (o'  -  a)  (R  sin  8  -  z,  -  Az ) 

bj  =  sec'  8*  (R  cos  8  -  x,  cos  G  -  Ax  cos  G  +  y,  sin  G  +  Ay  sin  G) 

G  =  o  -  y 

The  solutions  are  then  used  to  find  the  corrected  values  of  the  observer’s  coordi¬ 
nates. 


tan  (oc 
tan  4'  = 


Xj,  =  X,  +  Ax 
>0  =  y.  +  Ay 
Zg  =  z,  +  Az 

Procedurt 


Linear  Method 

The  given  information  consists  of  the  right  ascension  (o),  declination  (8),  and  the 
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radius  vector  (R)  for  each  observation  computed  from  orbital  data  together  with  the  right 
ascension  (o'),  declination  (8')  and  the  tirrie  of  each  observation  computed  from  photo* 
graphic  plates. 

Computation  of  y: 

The  hour  angle  (y)  of  the  Vernal  Equinox  is  computed  from  the  time  of  the  observa¬ 
tion  and  from  an  ephemeris  for  the  year.  A  description  of  the  method  is  included  in  the 
sample  computations  which  follow. 

Computation  of  x,,,  Yo,  and  z^: 

Tabulate  the  necessary  factors  to  be  used  in  the  equations.  Substitute  the  values 
into  the  equations,  writing  two  equations  for  each  observation.  The  linear  equations  can 
then  be  solved  for  x^,  yoi  and  Zo.  These  coordinates  may  then  be  converted  into  latitude 
and  longitude.  Appendix  H  contains  the  equations  for  this  conversion. 

Differential  Method 

The  given  information  is  the  same  as  that  for  the  Linear  Method  with  the  addition  of 
the  approximate  coordinates  (x,,  y, ,  z,)  of  the  observation  site. 

First  we  compute  and  8^  using  the  equations: 


tan  (a^  -  a) 


X,  sin  G  •»  y,  cos  G 
-X,  cos  G  +  y,  sin  G  R  cos  8 


.  cos  (a!  -  a)  (R  sin  8  -  z,) 

tan  8  - - ^ - r-7^ — ^  , 

-X,  cos  G  +  y,  sin  G  ->■  R  cos  8 

Then,  find  the  differences  between  these  computed  values  and  the  observed  values: 

Aa'  -  o'  -  a; 

A8'  =  8'  -  8i 

Find  the  eight  coefficients  for  the  differential  equations.  After  this  step  the  differ¬ 
ential  equations  can  be  formed  and  solved  for  the  increments  to  add  to  the  approximate 
coordinates. 

Usually,  several  approximations  will  be  necessary.  When  the  increments  are  added 
to  the  approximate  coordinates  of  the  first  computation,  the  resulting  values  can  be  used 
for  the  approximate  coordinates  for  the  second  approximation.  This  can  be  repeated  as 
many  times  as  necessary  until  the  increments  become  very  small. 

Sampit  Computation 

Computation  for  y  -  The  hour  angle  of  the  Vernal  Equinox: 

The  hour  angle  of  the  Vernal  Equinox  is  sidereal  time.  In  this  report,  our  reference 
point  is  the  Greenwich  meridian,  and  so  we  start  with  the  Universal  Time  (UT)  r»r  Green¬ 
wich  Civil  Time  (GCT)  of  the  observation. 

If  the  UT  of  observation  is  1%8  August  25  at  1  31*98,  we  must  find  the  side- 
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real  time  at  Greenwich.  From  The  American  Epbemeris  and  Nautical  Almanac  for  1958,  we 
find  that  Greenwich  sidereal  time  on  August  24  at  GCT  0  "  (midnight)  is  22  "11 00f488. 
The  solar  interval  of  1  SI*"  3l!98  must  be  converted  to  its  equivalent  in  sidereal  time. 
The  sidereal  day  is  S'"  56*556  longer  than  the  civil  day;  therefore  a  correction  must  be 
added  to  convert  solar  time  to  sidereal  time.  The  almanac  for  any  year  provides  tables  for 
easy  conversion. 

The  mean  solar  interval  of  I "  51*"  31*98  must  be  increased  by  18!322,  and  the  equiv¬ 
alent  sidereal  time  is  I**  SI*"  50!302. 

At  GCT  O"  the  sidereal  time  was  22"  11*"  00?488.  Hence  the  sidereal  time  of  the 
observation  is  24  "  02  "*50?79. 


Lineor  Method 


Given  Information: 

Observation  One 
a  =  17"  28'"  53!333 
S  =  40“  24'  09" 
R»  1.146813  e.r. 

14"  15'"  58!75 
«'=  39“  57'  08'.\)7 
y  -  24"  02'"  50!79 


Observation  Two 
a  =  6"  06'"  33* 

S  =  43“  32'  42" 

R  =  1.126957  e.r. 
a*  =  9"  06'"00!02 
8’  =  72“  38’  24'.'21 
y  -  11"  SB*"  09f47 


Tabulation: 


G' 

= 

14"  13"'  07?% 

G’ 

* 

2"  52"’  09.*45 

sin  G' 

= 

-0.548777237 

sin  G* 

= 

-0.682500802 

cos  G' 

= 

-0.835968627 

cos  G* 

= 

0.730884844 

sin  S 

= 

0.648153129 

sin  S 

0.688924071 

cos  S 

= 

0.761510027 

cos  8 

= 

0.724833515 

(o'  -  a) 

= 

-3"  12"’  54?583 

(a*  -  a) 

= 

2"  59"’27?02 

sin  (a*  -  a) 

= 

-0.745795003 

sin  (a*  -  a) 

= 

0.705408843 

cos  (a*  -  a) 

= 

0.666175513 

cos  (a*  -  a) 

= 

0.708800652 

cot  S' 

1.193772997 

cot  S' 

= 

0.312613361 
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Equations: 

-0.54P777237  Xo  -  0.835%8627  yo  =  -0.651311071 

-0.835%8627  Xo  +  0.548777237  -  1.193772997  Zo  =  -0.305566988 

-0.682500802  Xo  +  0.730884844  yo  -  0.576217590 

0.730884844  Xo  +  0.682500802  yo  -  0.312613361  Zo  -  0.336279009 

The  solution  of  these  equations  by  Doolittle’s  methcd  yields:** 
Xo  =  0.005833871  e.r. 

=  -37,209.63  meters 
yo  «  0  782937022  e.r. 

=  4993,733.93  meters 
Zo  =  0.619%9093  e.r. 

=  3954,290.84  meters 

The  geodetic  coordinates  of  this  point  are: 

^  =  38"  33'  45178 
A  =  90"  25'  36'.'91 

Difftrtntial  Method 


Given  Information: 

Observation  One 
o  =  17"  28"'  53?333 
8  =40°  24'  09" 

R  =  1.146815  e.r. 
a'-U*'  15"’  58!75 
8' =39"  57'  08.'07 
Y  =  24"  02"’  50?79 


X,  =  -0.005750000 

y,  =  0.789900000 

z,  =  0.619950000 


Observation  Two 
a  =  6"  06"*  33* 

0  =  43"  32'  42" 

R  =  1.126957  e.r. 
o' =9"  06"’00!02 
S'  »  70"  38'  24'.'21 
y  -  11"  58"’  09f47 
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Tabulation: 


G 

m 

17"  26'"  02?5i3 

G 

-5"  51"’  36M7 

sin  G 

m 

-0.989043181 

sin  G  > 

-0.999329646 

cos  G 

m 

-0.147626508 

cos  G  « 

0.036609552 

cos  8 

m 

0.761510027 

cos  8  « 

0.724833515 

sin  8 

m 

0,648153129 

sin  8  a 

0.688924071 

(o'  -  a) 

m 

-3"  12'"  54!583 

(o'  -  a)  - 

2"  59'"27?02 

sin  (c'  -  a) 

a 

-0.745795003 

sin  (o'  -  a)  = 

0.705408843 

cos  (a'  -  a) 

= 

0.666175513 

cos  (o'  -  o)  * 

0.708800652 

tan  (o'  -  a) 

* 

-1.119517287 

tan  (o'  -  o)  « 

0.99521472 

sec*  (a‘  -  a) 

= 

2.253318953 

sec*  (o' -  o)* 

1.99045235 

tan  8' 

as 

0.83768021 

tan  8'  > 

3.19883964 

sec*  8' 

s 

1.70170814 

sec*  s'  * 

11.23257504 

Computation  of  Aa*  and  A8'  for  Observation  One: 


tan  (ai  -  a) 


-0.110923180 

0.091217061 

-1.216035458 


oi  = 


cos (o’  -  a) 
tan  5c 


8'  = 


a*  -  aj.  *  Aa' 


8’  -  8i  -  A8' 


211®  39'  15'.'128 
0.635161721 
0.858991162 
40®  39'  44':234 
=  2®  20'  26'.'127 
=  0.040851016  radians 
--0?42'  36'.'164 
=  -0.012392632  radians 
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Computation  of  Aa'  and  AS'  for  Observation  Two: 


tan  (oj  -  o) 


0.034664031 

0.027696221 


-  1.251579789 

a;  =  143“  00'  46?755 
cos  (a*  -  a)  =  0.624213892 
tan  S,'  =  3.525775217 
Si  «  74®  09'  55'.'108 
a' -a;  -  Aa'  =-6“30'  46'.'455 

-  -0.113671622 

5  -  si  =  AS'  -  -1®  31'  30'.'898 

-  -0.26620623 


Coefficients  for  Observation  One: 


a,  -  -0.823772754 
aj  -  -1.254877447 

aj  =  0.205541131  i  0.332649558  Ax  -  2.228629851  Ay 

bi  =  -0.123663802 

bj  =  0.828501888 

bj  =  -0.666175513 

b4  =  -0.092002565  +  0.745795020  Az 

bs  =  0.155224808  +  0.251217183  Ax  -  1 .683062847  Ay 

Coefficients  for  Observation  Two: 


a,  =  -0.962895280 
aj=  1.031157127 

33  =  0.055127990  -  0.072869565  Ax  -  1.989118050  Ay 

b,  =  0.117108086 
bj-  3.1%695254 
bj  =  -0.708800652 

b*-  0.110352609  -  0.705408840  Az 

bj  -  0.311099780  -  0.411219518  Ax  -  11.225045283  Ay 
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Equations: 

-0.837361828  dx  -  1.163835657  dy  *  0.008396564 

-0.120550559  dx  t-  0.807644309  dy  -  0.6966419%  dz  -  -0.005682042 
-0.971178482  dx  ♦  0.805050853  dy  -  -0.006266490 

0.106161166  dx  +  2.897877564  dy  -  0.788985620  dz  -  -0.020825633 

A  solution  these  equations  produces  the  corrections  to  the  approximate  coordi* 
nates.  The  units  are  earth  radii  which  can  be  converted  to  meters  by  multiplying  by 
6.378.206.4. 

Ax-  0.000297675 
Ay  =  -0.007365136 
Az  -  -0.000536223 

Adding  these  algebraically  we  get: 

X,  ♦  Ax  =  -0.005452325 

y,  +  Ay  .  0.782534863 

z.  ♦  Az  -  0.619413777 

The  corresponding  geodetic  coordinates  for  this  first  approximation  are: 

0=38®  33’  0?.'695 
A  *9(f  23'  57.130 

These  values  were  then  used  for  x,,  y,  and  z,  in  a  second  approximation.  The 
resulting  corrections  and  coordinates  are: 


Ax  = 

-0.000379847 

Ay  = 

0.000402620 

Az  = 

0.000553518 

x,  +  Ax  = 

-0.005832172 

II 

+ 

■ 

0.782937484 

z,  +  Az  = 

0.619%7295 

0  = 

38^  33'  45'.'433 

A 

9(f  25'  36'.'457 
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A  third  approximation  results  ia  the  following  values: 

Ax  -  0.000000064 
Ay  .  -0.000001631 
Az  -  0.000001543 

X.  4^  Ax  «  -0.005832108 

y,  +  Ay  -  0.782935853 

z.  4^  Az  >  0.619968838 

^  -  38®  33'  4S'.'89 
A  -  90®  25'  36'.'44 

CONCLUSION 

The  preceding  materials  have  presented  the  mathematical  procedures  for  the  reduc¬ 
tion  of  artificial  satellite  observations.The  observation  of  these  satellites  offers  a  new 
opportunity  to  obtain  valuable  geodetic  information  not  previously  available,  and  will 
permit  the  unification  of  various  continental  datums  in  a  single  World  Geodetic  System. 

Observations  of  a  satellite  can  help  the  geodesist  to  determine  some  very  basic  and 
essential  data.  A  study  of  the  irregularity  of  the  orbit  can  reveal  data  concerning  the 
earth’s  gravitational  field.  From  investigations  of  a  satellite’s  period  and  acceleration, 
more  precise  data  can  be  collected  from  which  to  determine  the  size  of  the  earth.  The 
rate  of  change  of  the  plane  of  the  orbit  will  help  fix  the  flattening  of  the  earth.  Finally, 
the  satellite  provides  the  geodesist  with  a  direct  method  for  connection  datums. 

The  successful  use  of  earth  satellites  for  geodetic  purposes  requires  a  satellite  of 
high  stability.  It  also  requires  sufHcient  observation  stations  to  precisely  monitor,  record 
and  predict  the  satellite’s  orbit.  Auxiliary  observation  sites  can  then  be  used  to  gather 
the  desired  geodetic  data.  The  accuracies  obtainable,  as  shown  in  this  paper,  will  permit 
the  use  of  earth  satellite  observations  for  the  refinement  and  extension  of  the  present 
World  Geodetic  System. 
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APPENDIX  A 


% 


i 

r 


1 


Tti*  Ellips*  and  Sphorlcol  Right  Triongit 


1.  Ellipse 

a.  Axes 

OA  =  OA’  =  semimajor  axis 
denoted  by  the 
letter  a 

OB  »  OB'  =  semiminor  axis 
denoted  by  the 
letter  b 


Fig.  1 

b.  Eccentricity  e 

a2. 

e^  = - 2 — ;  The  value  of  e  determines  the  divergence  of  an  ellipse  from  a 

circle  and  varie '  from  0  to  1. 

That  is;  e  -0,  when  b>a; 

e  «  1,  when  b  ■  0. 

c.  Foci  F,  F' 

d.  cos  V 

a(l  -  e*)  -  R 

cos  V  - - 5 - 

e  n 

Let: 

KF-m,  FP-R 

Then, 

M  P  -  m  +  R  cos  B. 

Now,  as  given  by  Osgood  and 
Grausteii  ^ : 

Fig.  2 


B 


e  -  FP/MP 
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I 


W'V 


m  +  R  cos  9 


1  -  e  cos  9 

1, 

m  »  KO  -  FO,  where  by  Osgood  and  Graustein 
KO  «  a/e,  and  FO  -  C  -  ae 


Then 


m  -  —  -  ae  ■— (1  -  e*) 
e  e 


R  .  — — —• ,  but  6  -180  -  V,  or  cos  6  -  -  cos  v 

1  -  e  cos  9 


Therefore 


a  (1  -  e^ ) 

1  +  e  cos  V 


cos  V 


a  (1  -e*)  -R 


Spherical/Right  Triangle  and  Basic  Formulae 
a.  sin  a  -  sin  A  mn^ 


b.  tan  b  -  cos  A  tan  c 
w^ieref 

angle  C  is  a  right  angle 
and  a  »  BC,  b  -  AC,  c  »  ^ 


Fig.  3 
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APPENDIX  8 


Bosic  Eltments  of  a  Sefoilitt*s  Orbit 

1.  Kepler’s  Laws® 

a.  The  orbit  of  each  planet  is  an  ellipse  with  the  sun  at  one  of  its  foci. 

b.  The  radius  vector  of  each  planet  describes  equal  areas  (areal  velocity)  in 
equal  times. 

c.  The  squares  of  the  periods  of  the  planets  are  proportional  to  the  cubes  of 
their  mean  distances  from  the  sun. 

2.  Expression  for  the  areal  velocity 

Areal  velocity  is  defined  as  the  area 
swept  out  by  the  radius  vector  (r)  in 
unit  time.  By  Kepler’s  second  law 
above,  it  may  then  be  stated  that  the 
areal  velocity  of  an  earth  satellite  is 
constant.  Referring  to  Fig.  4,  this 
means; 

triangle  SFS'  *  triangle  S'  FS"  (in 
area),  if  the  time  interval  between  S 
and  S'  and  between  S'  and  S"  are 
equal.  F  (the  foci)  is  the  center  of  the 
earth.  S,  S' ,  S",  and  S'"  are  positions 
of  a  satellite  in  its  orbit.  The  r^ ,  r^  , 
r^ ,  and  r^  are  radius  vectors. 

In  Fig.  5: 

If  SS  ’  is  the  length  of  the  path  described 
in  unit  time,  then  SS'  is  the  linear  velocity  V. 

If  the  area  SFS’  is  the  areal  velocity  A,  then: 

A=  1/2  V  X 

where, 

X  =  perpendicular  from  F  to  the  line  of 
motion  (FB') 

V  *  SS'  or  the  base  of  the  triangle  SFS' 

We  next  find  expressions  for  V  and  X. 


Fig.  5 
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In  Fig.  6: 


The  earth  is  at  the  focus  F; 

FB  ”  A'  0  =  OA  -  a,  the  semimajor 
axis.  The  semiparameter,  P,  is  FY*  or 
FY.  The  eccentricity,  e,  is  represented 
by  the  expression  FO/A'  0.  From 
Russell,  Dugan,  and  Stewart’  we  have: 

V*-G(M+m)(---) 
r  a 

where, 

G  -  constant  of  gravitation  > 

6.668  X  10“*  dynes^ 

M  -  mass  of  the  earth  » 

5.983  X  10"'  gm.' 

m  -  mass  of  the  satellite 

r  -  radius  vector  measured  fron  F  to  any  point  on  the  orbit. 

When  the  satellite  is  at  the  extremity  of  the  minor  axis.  Point  B,  of  the  orbit 
(Fig.  6),  r  »  FB  >  a,  and  the  formula  for  V  becomes: 

V  yC  ni) 

From  Fig.  6,  it  is  seen  that  the  perpendicular  distance  X  for  the  satellite  at 
B  is  equal  to  FB' ,  but  this  is  evidently  equal  to  the  minor  axis  b  of  the  ellipse. 
Hence, 

A  =  1/2  VX  =  1/2  b 

Now,  by  the  geometry  of  the  ellipse:  b  *  a  Vl  -  e^  and  P  (the  semiparaineter)  = 

a(l-e2)  =  b-. 

a 

The  expression  for  the  areal  velocity  may  then  be  written: 

Vi 

A  *  1/2  P  VG  (M  +  m)  ,  that  is,  the  areal  velocity  is  proportional  to  the 

square  root  of  the  semiparameter,  and  since  the  areal  velocity  is  constant,  the 
above  formula  applies  to  any  position  of  the  satellite  in  its  orbit. 


Fig.  6 
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3.  The  expression  for  a  satellite’s  period 


The  whole  area  of  an  ellipse  is  «ab,  and  since  it  is  swept  nut  by  the  radius 
vector  at  the  rate  A,  the  period  will  be; 

p  _  a^^ 

~  A  "  VG  (M  +  m) 


let 


(M  +  m) 


Fig.  7 

4.  The  Eccentric  Anomaly  (E)® 

PS:0S  =  b:a  a  -  CA  .  CQ 

b  -  CD 

P  *  Position  of  satellite  in  its  elliptical  orbit 

Where: 


PS  -  R  sin  V,  QS  -  CQ  sin  E  »  a  sin  E 

R  -  FP  (radius  vector) 

R  sin  V  :  a  sin  E  B  b;a 

V  *  angle  AFP  (true  anomaly) 

or  R  sin  v  «  b  sin  E 

(1) 
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Now: 


FS  -  R  cos  V,  FS  «  CS  -  CF  ■  a  cos  F,  -  ae 

Then: 

R  cos  V  -  a  (cos  E  -  e)  (2) 

(1)*+  (2)*  ■  R^  sin*  V  -f  R*  cos*  v  «  b*  sin*  E  +  a*  (cos  E  -  e)* 

R*  ■  a*  (1  -  e*)  sin*  E  ♦  a  *(cos*  E  -  2  e  cos  E  +e*) 

R  *  -  a*  ( (1  -  e*)  sin  *  E  ♦  cos*  E  -  2  e  cos  E  +  e  *  1 

■  a*  (1  -  e  *  sin*  E  -  2  e  cos  E  +  e*) 

■  a*(l  +  2*(l  -  sin*  E)  -  2  e  cos  E 

■  a*(  1  +e*  cos*  E  -  2  e  cos  El 

R  *  a  (1  -  e  cos  E)  (radius  vector)  (3) 

Also: 


sin  1/2  V 


-  cos  V 


2  sin*  1/2  V.  1  -  cos  v 
2  R  sin*  1/2  v  *  R  (1  -  cos  v) 

»  a  (i  -  e  cos  E)  -  a  (cos  E  -  e) 

=  a  ( 1  -  e  cos  E  -  cos  E  +  e  1 

=  a  [(1  +  e)  -  cos  E  (1  +  e)  ] 

2  R  sin*  1/2  v  ■  a  (1  +  e)  (1  -  cos  E) 

2  R  cos*  1/2  V  =  a  (1  -  e)  (1  +  cos  E) 

by  division 


tan*  1/2 


1  +  e  1  -  cos  E 
1  -  e  1  +  cos  E 


see  (2) 


(similar  procedure) 


f  1  +  e  1  ^ 
,/2  -  l^j^J 


tan  E/2 


(4) 
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APPENDIX  C 

The  Earth's  Oblotenass  From  Solallitt  Obsorvotion 


The  mean  rate  of  rotation  of  the  orbital  plane,  or  the  mean  rate  of  change  in  the  right 
ascension  of  the  ascending  node,  0,  for  a  near  satellite  in  the  gravitational  field  is  repre¬ 
sented  by  the  following  theoretical  expression,  according  to  D.  G.  King-tlele;® 

nB^cos  i(J +J*BM||sinM  -1)  -j^DBM7sinM  .4)] 

where: 


n  SB 


the  mass  m  of  the  satellite  is  insignificant 


B 


a(l  -  e*) 


J  =  f-  —  -  f(---^^ 
J  2  '  ^2  r 


D.|(7f  -.5m) 

f  »  Oblateness  or  flattening  of  a  spheroid  representing  the  Earth 

m  =  the  ratio  of  centrifugal  to  gravitational  acceleration  at  the  equator,  multi¬ 
plied  by  (1  -  f).  The  value  of  m  is  accurately  known  to  be 
3449.79  X  10’*^  ±0.06  x  10"®. 


The  theoretical  expression  given  above  has  been  tested  against  observed  values  for 
the  change  in  Q  during  the  life-time  of  Sp  tnik  2.  It  was  found  that  the  observed  values 
were  always  slightly  lower  than  the  theoretical  with  the  ratio  between  them  remaining  vir¬ 
tually  constant.  The  value  for  the  flattening  used  was  1/297.1. 

Since  the  theoretical  values  involve  the  orbital  inclination  also,  tests  using  infor¬ 
mation  from  satellites  with  different  inclinations  are  necessary.  Preliminary  tests  with 
three  U.  S.  satellites  also  indicate  the  observed  values  of  the  change  in  Q  to  be  lower  than 
the  theoretical.  This  constant  discrepancy  seems  to  be  “due  mainly  to  the  J-term  and  that 
J  needs  to  be  reduced.’’® 

If  m  is  a  constant  then  the  J  function  depends  on  the  flattening  f.  Then,  if  J  is  to  be 
reduced  f  must  be  reduced.  In  the  case  of  Sputnik  2,  f  =  1/298.1  +0.1  would  have  been 
an  excellent  value. 
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APPENDIX  D 


Derivation  of  Formulae  for  Photographic 
Plate  Heductions 

In  order  to  develop  the  method  of  calculating  the  coordinates  of  an  unknown  point 
from  information  provided  by  the  photographic  image  of  a  portion  of  the  celestial  sphere, 
we  must  investigate  the  relationship  of  the  photographic  plate,  the  celestial  sphere  and 
the  plane  tangent  to  the  celestial  sphere. 


The  optical  axis,  CO,  is  perpendicular  to  the  photographic  plate.  The  length  CO  is 
the  focal  length  of  the  camera.  The  extension  of  the  optical  axis  intersects  the  celestial 
sphere  at  point  A  which  is  also  the  point  where  the  “Tangent  Plane”  is  tangent  to  the 
celestial  sphere.  The  tangent  plane  is  at  right  angles  to  CA  and  is  therefore  parallel  to 
the  photographic  plate.  Stars  at  L  and  have  tneir  projections  on  the  tangent  plane  at  N 
and  D  respectively.  The  image  of  L  appears  on  the  photographic  plate  at  M,  and  the 
image  of  B,  at  R. 
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By  angular  relationships,  we  can  see  that  the  system  of  images  on  the  photographic 
plate  is  similar  to  the  system  of  projections  on  the  tangent  plane.  The  difference  between 
the  two,  being  the  matter  of  linear  scale.  If  ^  =  OCM  ACL,  then 

^  OM  AN 
^  “6C  ’  AC 

Let  AlV,  AS'  be  the  positive  directions  of  the  rectangular  axes  in  the  tangent  plane, 
and  let  OR,  OS  be  the  positive  directions  of  these  axes  in  the  plane  of  the  plate.  OS  is 
parallel  to  AS'  and  OR  is  parallel  to  AR'.  Note  that  the  positive  direction  of  the  axes 
are  opposite  to  those  on  the  tangent  plane. 

If  we  let  i*  andij'  be  the  coordinates  of  the  projection  of  a  star  on  the  tangent 
plane,  and  i  and  i\  the  coordinates  of  the  image  on  the  plate;  then  by  similarity  we  have 


The  terms  i  and  ti  are  the  coordinates  of  the  image  on  the  plane.  In  order  to  deter¬ 
mine  the  relationship  between  the  C  and  tj  and  the  astronomical  coordinates  we  shall 
refer  to  Fig.  9  where  the  tangent  plane  and  celestial  sphere  are  tangent  at  A.  (The  camera 

points  towards  A.) 


/ 


axis 


T< - H' - >\ 


iCaxis 


Fig.  9 
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If  S  is  a  star  near  A,  its  projection  on  the  tangent  plane  is  obtained  by  extending 
the  line  CS  to  coincide  with  the  tangent  plane  at  T.  We  draw  the  great  circle  arc  AS 
which  projects  into  a  straight  line  AT  on  the  tangent  plane.  Let  P  be  the  north  pole  of 
the  celestial  sphere.  Then  AP  projects  into  AQ  which  we  shall  call  then'  axis.  We  draw 
AU  perpendicular  to  AQ,  and  let  this  be  the  axis.  The  positive  direction  of  ALi  is 
taken  to  be  eastward  from  the  meridian  AP,  so  that  increasing  values  of  correspond 
to  increasing  values  of  right  ascension. 

Since  the  arc  AP  projects  into  line  AQ  and  arc  AS  into  line  AT,  and  since  the 
sphere  is  tangent  to  the  plane  at  A  with  AT  X  to  CA,  then  QAT  defines  the  spherical 
angle  PAS. 

Let  arc  AS  be  denoted  by  ^  and  SAP  by  0;  then  QAT  ■  •.  Draw  perpendiculars  TU 
and  TV  to  AU  and  AQ  respectively.  Then 


VT-  ^'-  AT  sintf 

(4) 

UT  *  T|'  -  AT  cos  9 

(5) 

AT  -  AC  tan  ACT  -  AC  tan  6 

(6) 

£  ' 

*  tan  sin  9 

(7) 

=  tan  ^  cos  9 

(8) 

^  -  tan  ^  sin  9 

(9) 

-  tan  ^  cos  9 

(10) 

In  these  last  two  equations,  i  and  i{  are  the  coordinates  of  the  image  of  S  on  the 
photographic  plate  with  reference  to  rectangular  axes  thru  0,  the  center  of  the  plate,  and 
drawn  parallel  but  oppositely  directed  to  the  axes  AU  and  AQ  on  the  tangent  plane.  OC 
is  the  focal  length.  When  {  and  q  are  derived  and  in  the  same  units  as  the  focal  length, 
then  <f>  and  6  can  be  calculated.  Since  ^  and  0  are  functions  of  the  right  ascension  and 
declination  of  A  and  S,  the  right  ascension  and  declination  of  S  can  be  found  if  the  right 
ascension  and  declination  of  A  are  known. 

If  we  assume  the  unit  of  length  to  be  the  focal  length  and  (  and  tj  to  be  expressed 
in  this  unit,  then 

(  >  tan  ^  sin  $  (11) 

T) -tan  ^ cos 9  (12) 

(  and  T|  are  called  the  standard  coordinates  of  the  star  concerned.  The  use  of  these 
formulae  assumes  that  there  are  no  errors  involved  in  the  observations.  Since  a  perfect 
solution  is  only  approached,  never  reached,  a  discussion  of  the  errors  involved  and  their 
effect  on  the  formulae  will  be  included  at  the  end  of  the  discussion  of  the  derivation  of 
the  formulae. 
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Let  and  Do  be  the  right  ascension  and  declination  of  the  point  A  on  the  sphere, 
and  a  and  rt  the  coordinates  of  star  S.  The  relationship  between  i  and  n,  A*  and  D®, 
and  a  and  5  can  be  obtained  from  the  spherical  triangle  ASP  in  Fig.  9.  Arc  AP  =  90“  "  D, 
arc  SP  =  9(f  -  8,  ASP  =  a  -A,  arc  AS  =  SAP  =  6.  From  the  formulae  of  spherical 
trigonometry,  we  have 


cos^  =  sin  8  sin  CL  +  cos  8  cos  CL  cos  (a  -  Ao) 
sin  ^  sin  0  =  cos  8  sin  (a  -  Ao) 

sin  cos  0  =  sin  8  cos  CL  -  cos  8  sin  IX  cos  (o  -  Ao) 

Since  n  =  tan  ^cos  0  we  can  divide  (15)  by  (13)  to  get 

T|  =  cos  Do  -  cot  8  sin  Dp  cos  (a  -  Ap) 

sin  Do  +  cot  8  cos  Do  cos  (a  -  Ap) 

For  simplification  of  the  computations  let 
cot  q  =  cot  8  cos  (a  -  Ao) 

Then 

^  _  cos  Dp  -  sin  Dp  cot  q 
sin  Dp  +  cos  Do  cot  q 

or  11  =  tan  (q  -  CL ) 

Since  $  =  tan  ^sin  0  we  can  divide  14  by  13  to  get 

^  _ cot  8  sin  yg  -  Ap) _ 

sin  CL  +  cos  Do  cot  8  cos  (a  -  Ap) 


(13) 

(14) 

(15) 


(16) 


(17) 


(18) 


(19) 


using  q  cot  q  tan  (a  -  .Ap) 

^  sin  Dp  +  cos  Do  cot  q 

^  _  cos  q  tan  (a  -  Ap ) 
cos  (q  -  Dp ) 

derivation  of  (14) 


law  of  sines 


(20) 


sin  ^  cos  8 

sin  (a  -  A)  sin  0 


sin  ^  sin  0  =  cos  8  sin  (a  -  A) 
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derivation  of  (13) 

law  of  cosines  for  sides 

cos  if>  =cos  (90-8)  cos  (90-D)  +  sin  (90-8)  sin  (90-D)  cos  (a  -  A) 
cos  «l>  =sin  8  sin  0  +  cos  8  cos  D  cos  (a  -  A) 

derivation  of  (15) 

from  law  of  cosines  for  sides,  find  cos  6 

cos  (90-8)  =  cos  <ft  cos  (90-0)  +  sin  ^sin  (90-1))  cos® 
sin  8  -  cos  sin  0  =  sin  tf>  cos  D  cos  0 

sin  8  -  cos  if)  sin  D 
cos  0  =  :  T  TA 

sin  9  cos  U 

from  (14)  above 

,  CQS  8  sin  («  -  A) 

sin  <p  = - : — T - 

sin  0 

Multiply  sin  4>and  cos0 

,cos  8  sin  (g  -  A)^  ^inS  -  cos  <f)  sin  D 


sin  if)  COS0  ={ 


sin  COS0 


sin  0  ^  ^  sin  cos  0 

cos  8  sin  (a  -  A)  (sin  8  -  cos  sin  0] 


1 


sin  0  sin  ^  cos  D 

since  cos  8  sin  («  -  A)  is  equal  to  sin  0  sin<^,  by  cancelling,  we  have 

-  sin  8  -  cos  6  sin  0 

sin  d>  COS0  = - =r - 

cos  U 

substituting  the  value  of  cos  9 

sin  8  -  sin  D  [sin  8  sin  D  cos  8  cos  D  cos  (a  -  A)) 


sin  4)  COS0  = 


cos  0 


sin  8  -  sin  8  sin^  D  -  sin  D  cos  8  cos  D  cos  (a  -  A) 


cos 


D 


sin  8(1-  sin^  Dl  -  sin  D  cos  8  cos  D  cos  (a  -  A) 


cos  D 


sin  ^  COS0  =  sin  8  cos  D  -  sin  D  cos  8  cos  (a  -  A) 
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To  get  the  formulae  in  terms  of«,  8,  Aq,  Do,  i  and  tj  from  (16) 

cos  Do  -  cot  8  sin  Dp  cos  (a  -  Ap) 
sin  Dp  +  cot  8  cos  Dp  cos  (a  -  Ap) 


T)  sin  Do  +  11  cot  8  cos  Dp  cos  (a  -  Ap)  =  cos  Dt,  -  cot  8  sin  l\  cos  (a  -  Ap) 

cot  8  cos  (a  -  Ap)  [  T)  cos  Db  +  sin  Dip  1  *  cos  Do  -  tj  sin  EX) 

cos  Do  -  11  sin  Do 


cot  8  cos  (a  -  Ap)  = 
cot  8  cos  (a  -  Ap)  = 


from  (10) 


i 


TJ  cos  (X)  +  sin  Do 

1  -  TJ  tan  Dp 
TJ  +  tan  Dp 

cot  8  sin  (a  -  Ap) 


(21) 


sin  IX>  +  cos  Dp  cot  8  cos  (a  -  Ap) 

i  [sin  IX,  +  cos  Dp  cot  8  cos  (a  -  Ap)l  =  cot  8  sin  (a  -  Ap) 

.  .  /  *  V  .  r  .  cos  Do  (1  -  TJ  tan  Dp), 

cot  8  sin  u»  -  Ap)  =  t  [ sin  Do  + - - 1 

TJ  +  tan  l\> 

sin  rU 

„  r  •  n  r»  ^  cos  Do  COS  , 

=  t  f  sin  IX>  +  COS  Dp - p; - J 

TJ  tan  IX 

,  -  .  p,  cos  Db  -  sin  Dp  , 

[sin  Do  + - - — - 1 

TJ  +  tan  Dp 

TjsinDp  +  sinDbtanDp  +  cosDb-TjsinfX. 


(  [ 


TJ  +  tan  Do 


1 


sin^  fX  sec  Dp  +  cos  IX 

I - ^^ti^iTo; - ^ 

^  sin  Dp  tan  Dp  +  cos  Dp 

I - ~lanTX - ^ 


.  8  •  r  AT  4  sec  Dp 

cot  8  sin  (a  -  Ap)  - - 

TJ  4  tanDp 


(22) 
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Divide  (22)  by  (21) 


cof  S  sin  (d  -  Ao) 


cot  8  cos  ip  -  \q) 


i  sec  l\ 

T|  tan  Db 
1  -  T|  tan  Db 
1)  4  tan  Db 


tan  (a  -  A®) 


i  sec  Db 
1  -  1)  tan  Do 
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APPENDIX  E 


Discussion  of  Errors  Portoining  to  tho  Photogrophic  Method 

The  standard  coordinates  differ  from  the  plate  coordinates  by  small  amounts  which 
represent  measurable  or  calculable  amounts  of  imperfection  due  to  a  number  of  errors  in* 
volved.  These  errors  fall  into  two  categories: 

1.  Geometrical  and  Mechanical 

(a)  Error  of  orientation 

(b)  Non*perpendicularity  of  axes 

(c)  Centering  error 

(d)  Error  of  tilt 

2.  Spherical 

(a)  Refraction 

(b)  Aberration 

Taken  separately,  these  errors  produce  in  each  coordinate  a  displacement  of  the 
image  of  a  star  on  the  plate  from  the  position  corresponding  to  the  standard  coordinates 
of  the  star.  This  displacement  can  be  represented  as  a  linear  expression  since  a  study  of 
these  errors  reveals  that,  generally  second  order  effects  are  negligible. 

The  linear  expression  obtained  from  the  combinatio.i  of  all  the  effects  on  the  coordinates 
is: 

^  -x  =  ax  +  by  +  c  (24) 

ij  -  y-dx+ey+f  (25) 

The  terms  a,  b,  ...,  f  are  called  plate  constants.  The  values  reflect  the  errors  pre¬ 
sent  in  any  one  plate. 

Since  the  x  and  y  coordinates  differ  from  the  (  and  rj  coordinates  by  small  amounts, 
the  plate  constants  can  be  computed  from  the  plate  coordinates  or  from  the  standard  co¬ 
ordinates.  The  plate  coordinates,  coming  from  observational  data,  have  more  avenues  for 
errors  than  the  computed  standard  coordinates;  therefore,  it  is  recommended  that  the  plate 
constants  be  computed  from  the  form: 

f  -x-a^+bii+c  (26) 

V  -  y=d(+ei}+f  (27) 
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APPENDIX  F 


Development  of  the  Method  of  Dependences 

In  the  working  procedure,  the  method  of  solving  for  the  plate  constants  by  the  least 
squares  is  explained.  The  process  requires  quite  a  lengthy  computation,  and  if  we  have 
several  plates,  the  plate  constants  have  to  be  determined  from  each  plate  which  requires 
a  large  amount  of  numerical  work.  Consider  a  series  of  plates  with  the  same  plate-center 
and  in  which  the  position  of  the  unknown  point  differs  little  from  one  plate  to  the  next. 
The  method  of  dependences  was  developed  mainly  for  use  with  such  a  series  of  plates. 

It  works,  of  course,  with  one  plate  but  nothing  is  gained  in  the  way  of  shortening  the 
computation  unless  it  is  applied  to  several  plates. 

The  method  of  dependences  eliminates  the  specific  values  of  the  plate  constants 
and  replaces  them  with  terms  which  express  their  functions. 

The  solution  will  be  the  same  as  given  in  the  working  instructions  up  to  the  point 


of  solution  of  the  normal  equations  which  are: 

a  +  bfi,,,.]  +  c  ii,]  =  -x,)l  (1) 

a  +  bfijjiij]  +  c  [ijjl  =  -Xj)l  (2) 

a  [^jl  +  bijfj]  +  nc  =  (3) 

For  the  unknown  point  we  have  t  •  x  *  a{  +  biy  +  c  (4) 


If  (^t  »  standard  coordinates  of  the  unknown  point  for  one  of  the  plates 

(which  we  shall  refer  to  as  the  “selected  plate”),  then  we  can  write  (4)  as  follows: 

{  -  X  =  a{,  +  bij,  +  c  +  a(^  -  4,)  +  Utj  -  ij,)  (5) 

It  has  already  been  pointed  out  that  the  quantities  ({  -  {,)  and  (tj  -  tj,)  are  very 
small  quantities  and  the  values  of  a  and  b,  are  al.  ^  small.  Therefore,  we  can  neglect  the 
last  two  terms,  and  the  equation  for  the  unknown  point  becomes: 

{  -  X  =  a^j  +  bij,  +  c  (6) 

We  can  manage  to  eliminate  a,  b,  and  c  from  the  four  equations  (1),  (2),  (3)  and  (6). 


Multiply  these  in  order  by  P,  Q,  R,  and  -1. 

Then  if  P,  Q  and  R  are  given  by: 

P  U.ijJ  UJ- 

(7) 

P  »»(!  +Q  [>?,»?,]  +  R  >», 

(8) 

P  U,1  +  Q[l,l  +Rn  “  1 

(9) 
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we  have: 


^  -  X  -  -  x,)l  +  0tili({i  -x,)]  +  R[^,  -X,  J 

which  can  be  written 
n 

{-X  -2  (P^.+OiJ.+RXe,  -x.)  (10) 

i  -1 

let  D,  -P^, +0i|, +  R  (11) 

then  n 

{  -X  -  2  D,  (f,  -X  )  (12) 

i  -1 


Since  x  differs  from  ^  by  a  small  amount  we  can  substitute  the  measured  coordinates 
for  the  standard  coordinates.  Thus  P,  0  and  R  will  be  determined  from  (13),  (14)  and 
(15). 


P[x,x,]  +  0lx,y,  1  +  R[x,]  *Xj 

(13) 

PI  x,y,l  +  0l  yiyjl  +  Rl  y,l  -  y, 

(14) 

Pi  x,  ]  +  Q 1  y,  ]  +  Rn  =  1 

(15) 

D,  =  Px(  +  Qyj  +  R 

(16) 

The  dependences  can  be  formed  from  (16)  and  checked  by: 


2  D,  =  I  (17) 

i  =  1 


Then  the  solution  for  the  standard  coordinate  {  can  be  found  from  (12). 
The  procedure  is  tlie  same  foriy. 
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APPENDIX  G 


Development  of  Formulae  for  the  Geocentric  Coordinates 
of  the  Observation  Station 

In  order  to  derive  the  formulae  which  define  the  relationship  between  an  observer’s 
position  and  the  position  of  a  small  earth  satellite,  we  must  relate  each  to  an  eartli- 
centered  set  of  axes. 

Let  the  axes  X,  Y,  Z  be  defined  so  that  the  Z  axis  points  toward  the  north  pole, 
the  X  axis  points  toward  the  vernal  equinox  and  the  Y  axis  is  six  hours  eastward 
from  the  vernal  equinox.*® 

Another  set  of  axes,  x,  y,  z  are  centered  at  the  same  point,  and  the  z  axis  coincides 
with  the  Z  axis  of  the  first  system.  The  x  axis  points  toward  the  Greenwich  meridian, 
and  the  y  axis  is  located  90°  westward  from  the  Greenwich  meridian.** 

Consider  the  observer’s  position,  K,  at  a  distance,  p,  from  the  center  or  origin,  C. 

In  order  to  obtain  the  geodetic  coordinates  x«,  y^ ,  z^,  we  get  the  relationships  from  the 
geometry.  Let  K'  be  the  projection  of  K  in  the  xy  plane.  Then  the  distance  that  CK '  is 
from  the  x  axis  is  the  angle  A* .  Then  let  be  the  angle  between  CK*  and  CK.  In  other 
words  ,4>' ,  A' ,  and  p  are  the  geocentric  coordinates  of  K  in  polar  form. 

Then  Zo  >>  p  sin  ^ ' 

yo  =  p  cos^'  sin  A* 

Xj,  -  p  cos^'  cos  A'  (1) 

If  Xj ,  Yj ,  Z I  are  the  coordinates  of  K  in  the  celestial  system,  we  find 

Z|  >  p  sin 

Yj  -  p  cos^'  sin  6 

X  j  «  p  cos^i  cos  9  ^2^ 

where  9  is  the  hour  angle  of  the  vernal  equinox  with  respect  to  the  geocentric  meridian 
at  K. 

The  hour  angle  of  the  vernal  equinox  with  respect  to  the  geocentric  meridian  at 
Greenwich  is  y. 

Then  y  -  ^  ♦  A' 
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Substituting  in  (2)  we  get 


'f, 

= 

p  cos 

cos 

(y-X‘) 

p  cos 

(cos 

;  y  cos  X* 

4  sin  y  sin  X' ) 

- 

Xo  cos  y 

♦  Yo 

sin  y 

Y 

1 

ar 

p  cos 

sin 

(y  -  X') 

Y, 

JS 

p  cos 

(sin 

y  COS  X* 

-  cos  y  sin  X') 

Y. 

Xo  sin  y 

-Yo 

COS  y 

Z,  =  z„  (3) 

If  a  and  £  are  the  geocentric  rif^t  ascension  and  declination  of  the  satellite,  and 
R  the  distance  from  the  center  of  the  earth  to  the  center  of  the  satellite,  then  the  celestial 
geocentric  coordinates  of  the  satellite  are 

*  R  cos  S  cos  a 

Y  *  R  cos  S  sin  a 
2 

Zj  =  R  sin  6  (4) 

If  r  is  the  distance  from  K  to  S  and  o'  ,  are  the  apparent  coordinates  of  S  as 
seen  from  K,  we  have  from  geometric  relationships 

Xj  -  ^ '  cos  a ' 

Y.  -  Y,  =  r  cos  S'  sin  o' 

t  1 

-  Z|  =  r  sin  S'  (5) 

These  equations  are  the  basic  ones  which  relate  the  apparent  coordinates  of  the 
satellite,  the  geocentric  coordinates  of  the  satellite,  and  the  geocentric  coordinates 
of  the  observer’s  position. 

We  must  put  equations  (5)  into  a  form  which  uses  the  observed  data  and  computed 
data  and  can  be  solved  for  the  observer’s  coonlinates. 

If  we  multipl,  i  sides  of  the  first  equation  of  (5)  by  (  +  sin  a '  )  we  have 
sin  «'  (Xj  -  Xj )  =  sin  a'  r  cos  6'  cos  a' 

Now  if  we  substitute  the  values  for  X^  and  Xj  from  (4)  and  (3)  respectively  we  get 

sin  [  (R  cos  fi  cos  a)  -  (x^  cos  y  *  Vo  sin  y)  1  =  sin  a  '  r  cos  S'  cos  o' 
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Next  we -lake  the  second  equation  of  (5)  and  multiply  both  sides  by  (•  cos  a') 

-  cos  o*  (Yj  -  Yi)  *  -  cos  a*  r  cos  6'  sin  o* 
and  substituting  from  (4)  and  (3) 

-  cos  o'  [  (R  cos  5  sin  a)  -  (xosin  y  -  yoCOs  y )  ]  =  -  cos  o'  r  cos  6'  sin  o' 

We  multiply  both  sides  of  the  third  equation  of  (5)  by  sin  o'. 

Combining  these  two  equations  algebraically  we  have 

sin  o*  [  (R  cos  S  cos  o)  -  (XoCosy+  yo  sin  y)  ]  -  cos  o'  [(R  cos  S  sin  a) 

-  (Xgsin  y  -  yo  cos  y )  ]  =  sin  o'  r  cos  6'  cos  o'  -  cos  o'  r  cos  8'  sin  o* 

sin  o'  R  cos  8  cos  a  -  sin  o'  XgCos  y  -  sin  o'  yo  sin  y  *  cos  o'  R  cos  8  sin  o 

+  cos  o'  Xosiny  -  cos  a'  yocosy  =  sin  a*  r  cos  8  cos  a*  -  cos  o'  r  cos  8 'sin  o' 

Combining  terms  we  get 

Xosin  (a'-y)  +  yoCos  (a*  -  y)  =  R  cos  8  sin  (a*  -a)  (6) 

Again  using  the  equations  in  (3),  multiply  the  first  equation  by  (  +  cos  o'),  the 
second  by  ( +  sin  a* ),  and  the  third  by  ( -  cos  8* ) 

cos  a'  r  cos  8'  cos  a*  =  cos  a*  R  cos  8  cos  a  -  cos  o'  Xo  cos  y  -  cos  o'  yo  sin  y 

sin  o'  r  cos  8'  sin  a*  ~  sin  o'  R  cos  8  sin  a  >  sin  o'  Xq  sin  y  4'  sin  a*  yo  cos  y 

-  cot  8'  r  sin  8'  =  -  cot  8'  R  sin  8  +  cot  8'  Zo 

r  cos  8'  (sin^  o'  +  cos^  a')  -  r  sin  8'  cot  8'  =  R  cos  8  (sin  a*  sin  a  +  cos  a*  cos  a) 

•  R  sin  8  col  8'  -  Xo(sin  o'  sin  y  +  cos  o'  cos  y)  +  yo(sin  o'  cos  y-  cos  o'  sin  y)  +  ZoCot8' 

r  cos  S'  -  r  sin  S'  cos  8*  “  R  cos  8  cos  (o'  -  a)  -  R  sin  8  cos  8'  -  x^cos  (a*  -y) 

+  yoSin  (rf  -y)  +  Zo  cot  8* 

R  sin  8  cot  8'  -  R  cos  8  cos  (o'  -a)  ~  XoCos  (a*  -  y)  +  y©  sin  (o'  -  y)  +  Zo  cos  8' 

Or 

Xo  COS  (o'  -  y)  -  yo  sin  (o'  -  y)  -  Zo  cos  8'  *  R  cos  8  cos  (a'  -  a)  -  R  sin  8  cot  8'  (7) 

If  we  let  G'  *  o'  -  y,  then  equations  (6)  and  (7)  become 
Xo  sin  G'  +  yo  cos  G'  =  R  cos  8  sin  (o'  -  a) 

XoCos  G'  -  yo  sin  G'  -  ZoCot  8'  *  R  cos  8  cos  (o'  -  a)  -  R  sin  8  cot  8' 
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These  equations  will  produce  the  x,  y,  z  coordinates  of  a  station  which  have  before 
been  unknown.  In  the  case  of  a  station  that  is  located  by  approximate  coordinates,  it  is 
possible  to  provide  a  solution  which  gives  corrections  to  the  position.  The  development 
of  the  formulae  begins  again  with  the  equations  from  (5).  First  we  multiply  each  respect¬ 
ively  by  (-  sin  a),  (+  cos  a)  and  (0). 

•  r  cos  8*  cos  a*  sin  e  »  -  R  cos  S  cos  a  sin  a  +  XnCOs  y  sin  a  +  yosin  y  sin  a 
r  cos  s'  sin  a*  cos  a  -  R  cos  5  sin  a  cos  a  -  Xosin  y  cos  a  +  y©  cos  y  cos  a 
Combining  these  equations,  we  have 

r  cos  s'  sin  a  cos  a  •  r  cos  S'  cos  a*  sin  a  »  R  cos  S  sin  a  cos  o 
-  R  cos  6  cos  a  sin  a  +  XoCos  y  sin  a  -  x®  sin  y  cos  a  +  y^  sin  y  sin  a  +  y©  cos  y  cosa 

r  cos  S'  (sin  a'  cos  «  -  cos  a'  sin  a)  -  Xo(cos  y  sin  a  -  sin  y  cos  a) 

+  yo(  sin  y  sin  a  +  cos  y  cosa) 
r  cos  s'  sin  (a'  -  a)  *  x©  sin  (o  -  y)  +  y©  cos  (a  -  y)  (8) 

Next  we  multiply  each  equation  from  (S)  by  (-t-cos  a  ),  (+  sin  a)  and  (0) 
r  cos  s'  cos  a'  cos  a  =  R  cos  S  cos*  a  -  x©cos  y  cos  a  -  y©  sin  y  cos  a 
r  cos  S'  sin  a*  sin  a  >  R  cos  S  sin*  a-  x©sin  y  sin  a  +  y©cos  y  sin  a 
Combining  we  get 

r  cos  S*  (cos  tt'  cos  a  +  sin  o'  sin  a)  ■=  R  cos  S  (cos*  a  +sin*  a) 

-  x©(cos  y  cos  a  +  sin  y  sin  a)  +  y©(cos  y  sin  a  -  sin  y  cos  a) 


Or 


r  cos  S'  cos  (o'  -  a)  =  R  cos  S  -  x©cos  (a  -  y)  +  y©  sin  (a  -  y)  (9) 

Dividing  (8)  by  (9)  and  clearing  fractions 

r  cos  S'  sin  (  tt'  -  a)  x©  sin  (a  -  y)  -f  y^  cos  (a - y) _ 

r  cos  s'  cos  (tt'  -  tt)  X©  cos  (a  -  y)  +  y©  sin  (a  -  y)  +  R  cosS 

tan(«'  -  «)  ^ 

-X©  cos  (tt  -  y)+  y©  sin  (o-y)  +  n  cos  S 

-  x©cos  (tt  -  y)  tan  (a'  -  a)  +  y©sin  (a  -  y)  tan  (a'  -  a)  +  R  cos  S  tan  (a*  -  a) 

-Xosin  (a-y)  +  y©cos  (a-y) 

x©  [  cos  (tt  -  y)  tan  (  tt'  -  «  )  +  sin  (a-y)  1 

+  y©  [  -  sin  (tt  -  y)  tan  (  tt'  -  a  )  +  cos  (o-  y)l  ■  R  cos  8  tan  (a'  -  a)  (10) 
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Dividing  tho  third  equation  of  (5)  by  (9) 


r  sin  5'  cos  (a*  -  a  )  (zg  -  Zi) _ 

r  cos  T  R  cos  8  -  Xo  cos  (a  -  y)  +  sin  (a  -  y) 

cos  (a*  -  a  )  (l{  sin8  >  z.) 

tanS'  -  - - - r - ^-^-2 — 5 — - 5 

-  X,,  cos  (a  -  y)  +  yo  sin  (a  -  y)  +  If  cos  a 

-  cos  (o-  y)  tan  5'  +  y^  sin  (a  -  y)  tan  8'  +  H  cos  8  tan  8' 

-  cos  (a*  -  a )  (R  sin  8  •  z^) 

Xg  cos  (a  -  y)  tan  8*  -  sin  (a  .  y)  tan  8'  •  R  cos  8  tan  8' 

■  -  cos  (o'  -  o  )  (R  sin  8  -  Zg)  (11) 

If  we  let  G  -  a-y,  then  equations  (10)  and  (11)  become: 

Xg  I  sin  G  -«•  tan  (  a'  •  a )  cos  G  1  +  yg  I  cos  G  -  tan  (  a*  •  a )  sin  G  ] 

-R  cos  8  tan  (a*  -  a)»0  (12) 

x^  cos  G  tan  8'  -  yg  sin  G  tan 8'  -  Zg  cos  (a'  -  a) 

+  Rcos(a'  -  o )  sin  8  -  R  cos  8  tan  8*  *0  (13) 

I.et  X,,  y^  ,  z,  be  provisionally  adopted  coordinates  of  K  so  that 
Xg  -  X,  +  Ax,  yg  -  y,  +  Ay,  Zg  -  z,  +  Az 

Let  ,  8g  be  the  apparent  coordinates  obtained  through  use  of  the  values  x,, 
y,  ;  z,  above  so  that 


a'  -  Oc  +  Aa'  ^ 

8' 

•  8i^  +  A8 

or 

Aa'  -  a'  -  a^  , 

A8' 

-  8'  -  8' 

In  the  following  development,  we  will  assume  the  differentials  dx,  dy,  dz,  da'  and 
d6'  to  equal  Ax,  Ay,  Az,  Aa'  and  A8'  respectively,  giving: 

Xg  -  X,  +  dx,  yg  -  y,  +  dy,  Zg-  z,  +  dz 

do*- o' -  aj,  dfi'-6' .  fig'  (14) 

Differentiating  equations  (12)  and  (13)  with  respect  to  Xg,  yg ,  z^,  a'  and  fi', 
which  may  vary,  we  obtain; 
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a,  =  sin  G  +  tan  (o'  -  a)  cos  G 
32  =  cos  G  -  tan  (n'  -  a)  sin  G 

a^  =  sec^  (n'  -  n)  (  R  cos  S  -(x,  +  dx)  cos  G  +  (y,  +  dy)  sin  G  1 

bj  *  tan  8'  cos  G 

bj  =  -  tan  8'  sin  G 

bj  -  -  cos  ( n'  -  n  ) 

b^»sin(n'  -  n  )  I  R  sin  8  -  (z,  +  dz)  J 

b  *  sec^  8'  (  R  cos  8  -  (x,  +  (ix)  cos  G  +  (y,  +  dy)  sin  G  ] 
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APPENDIX  H 


Relationship  Between  Geocentric, 

Geodetic,  and  x,y,z  Space  Coordinate^ 

I.  The  X,  y,  z  space  coordinate  formulae  are: 

X  -  p  cos  ^ '  cos  A 
y  -  p  cos  ^ '  sin  A 
z  -  p  sin 

where: 

p  -  a  (l-e*  sin*  fl)** 
tan  “  (l"e*)  tan  ^ 

II.  Explanation  of  Notations 

1.  A  -  geodetic  coordinates  (latitude  and  longitude) 

2.  p  >  geocentric  radius 

3.  •  geocentric  latitude 

4.  a  •  semimajor  axis  of  the  reference  ellipsoid 

5.  e  -  eccentricity  of  the  reference  ellipsoid 

6.  0  *  reduced  latitude 

tan  0  -  (1-f)  tan  ^ 

7.  f  -  flattening  of  the  reference  ellipsoid 

8.  The  xz  plane  contains  the  Greenwich  meridian 

The  xy  plane  is  the  equatorial  plane;  the  positive  y-axis  being  90®  west 
of  the  x*axis. 

III.  Elements  of  the  Clarke  (1866)  ellipsoid  for  use  with  stations  on  North  American 
Datum 

a  -  6,378,206.4  meters 
e*  -  0.00676866 
f  "  0.00339008 
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IV.  To  find  geodetic  coordinates  from  x,  y,  z  space  coordinates 


tan  k  «  y/x 


tan^ 


tan  si* 
l-e^ 


where: 

sin  si*  *  z/p 

p  -  (x^  +  y2  +  z^)** 

V.  Numerical  Example 

Given: 

SS  -  35®  N 
A  -  80®  W 

Assume  the  point  to  be  on  North 
American  Datum,  then: 

a  *  6,378,206.4  meters 
e2  -  0.00676866 
f  -  0.00339008 

Tabulation: 

tan  -  0.70020755 
1  -  e^  -  0.99323134 
tan  -  0.69546808 
sin  -  0.57096264 
cos  si'  »  0.82097604 
sin  A  -  0.98480775 
cos  A  -  0.17364818 

Computation  of  x,  y,  z 


♦  Z 


Fig.  10 

1  -  f  -  O.9O06O992 
tan  «  «  0  .39783379 
sin  e  -  0.32749238 
sin^  e  -  0.32749238 
e2  sin2  e  -  0.00221668 
(l-e2  sin2  $)’*  -0.998891045 

p  -6.37113326  X  10* 


3,637,679.07  meters 


X  -  5.23054775  x  10*x0.17364818  -  908,275.10  meters 
y  -  5.23054775  x  10®  x 0.98480775  -  5,151,083.96  meters 
z  -  6.37113326  x  10®x  0.57096264 


VI.  Reverse  Check 


Given: 


X  »  908,275.10  meters 
y  -  5,151,083.%  meters 
z  •  3,637,679.07  meters 


Computation  of  A,  p,  and 

tan  A  -  y/x  -  5.07128171 
A  -80®  W 

sin  -  z/p  -  0.570%264 
tan  .  0.69546809 
l/(l-e2)  -  1.00681479 
tan  ^  -  0.70020756 
^  -  35®  N 


♦  y2  ♦  z2 

-  0.4059133864  x  10** 
p  -  6.37113323x  10* 
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